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-^superseded  by  later  SCOUR  calculations,  provide  an  excellent  in- 
troduction to  the  predictive  capacity  of  the  SCOUR  code. 

A comprehensive  comparison  is  made  between  the  SCOUR  code  pre- 
dictions and  experimental  observations  of  the  near-surface  dust 
clouds  produced  by  the  MIDDLE  GUST  II  (100  tons  TNT)  and  MIXED 
COMPANY  (500  tons  TNT)  high  explosive  detonations.  Included  are 
the  results  of  the  most  recent  (1975)  SCOUR  calculations  that 
were  made  with  an  improved  version  of  the  SCOUR  code.  These  re- 
sults indicated  that  the  nominal  SCOUR  prescription  of  dust 
lofting  and  diffusion  parameters  appear  to  be  appropriate  for 
the  positive  and  early  negative  phase  of  these  HE  blast  wave 
events.  At  later  times,  the  particle  impact/deposition  model 
in  SCOUK^lays  a major  role  in  the  level  of  predicted  concentra- 
tions. BiNad justing  model  parameters,  it  was  possible  to  explain 
the  discrep^cies  between  SCOUR  results  and  the  5 mm  glass  beads 
particulate  lofting  experiment  (MIDDLE  GUST  II) . This  suggests 
that  impact  and  deposition  phenomenology  should  be  further 
studied. 

The  last  area  of  investigation  was  the  application  of  a new  S^ 
code,  FLUB,  to  the  calculation  of  the  phenomenon  known  as  "re- 
verse percolation".  This  is  the  condition  that  occurs  when  the 
air  in  porous  soils  is  trapped  and  pressurized  during  the  early 
positive  phase  of  the  blast  wave  flow  field.  The  trapped  'air  is 
capable  of  rapidly  percolating  up  through  the  soil  matrix  when 
the  overpressure  is  reduced.  The  present  results  demonstrate 
that  with  a one-dimensional  treatment  of  the  problem,  porous 
(strengthless)  soils  can  be  fluidized  down  to  a depth  of  a meter, 
and  significant  quantities  of  sand/dust  will  be  raised  above  the 
original  surface.  It  is  speculated  that  this  can  affect  the  mass 
lofting  of  the  blast  wave  flow  field,  but  due  to  more  pressing 
research  and  development  problems  on  precursor  dust  storms,  the 
study  was  terminated. 
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I.  INTRODUCTION  AND  SUMMARY 


Systems,  Science  and  Software  (S^)  has  developed  a 
computer  code  (SC0UR2)  that  calculates  the  time  varying 
structure  and  intensity  of  the  near-surface  dust  storm 
produced  in  the  aftermath  of  nuclear  detonations.  This 
information  is  required  by  systems  designers  concerned  with 
"hardening"  missile  sites  to  withstand  the  high  dust  dynamic 
pressures  in  the  positive  phase  of  the  blast  wave  flow  field, 
and  for  those  concerned  with  near-surface  dust  concentrations 
in  the  late  time  (negative  phase)  dust  environment. 

SC0UR2  was  developed  over  a five  year  period  from 
December  1970  to  December  1975.  Most  of  the  effort  has 
been  reported  in  References  1-4.  This  report  completes  the 
documentation  of  SC0UR2  development  and  calculations.  It 
covers,  primarily,  the  third  annual  report  period  from  March 
1973  to  February  1974,  but  also  includes  some  unreported  work 
from  later  periods. 

The  SCOUR  code  is  briefly  described  in  the  second 
chapter  of  this  report.  An  appendix  is  provided  that  lists 
code  parameters,  nomenclature,  and  input-output  formats. 

The  latter  are  taken  from  the  most  recent  version  of  the 
code. 

Our  primary  interest  in  this  time  period  was  the  pre- 
diction of  the  nuclear  dust  storm  cloud  dimensions,  concen- 
trations, and  dynamic  pressure.  Approximately  fifteen  cal- 
culations were  performed  to  provide  state  of  the  art  input 

data  for  particular  detonations  that  were  of  interest  to  the 

[5-7] 

defense  community.  The  results  of  these  calculations 

are  summarized  in  the  third  chapter  of  this  report.  They 
include  surface  burst  calculations  from  32.5  KT  to  20  MT,  and 
the  results  of  the  first  SCOUR  calculation  with  a precursor 
flow  field  (36.6  KT  at  350  ft  HOB). 
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It  should  be  noted,  however,  that  the  results  pre- 
sented in  Chapter  III  have  been  superseded  by  the  more  re- 
cent calculations  reported  in  References  3-4.  The  latter 
were  necessary  because  a number  of  minor  coding  errors  in 
the  original  SCOUR  computer  program  were  detected  towards 

[31 

the  end  of  the  third  report  period.  Nevertheless,  these 

earlier  results  (reported  in  Chapter  III)  are  very  close  to 

[3  41 

those  calculated  with  the  SC0UR2  code  ' and  they  were 
crucial  to  the  continued  development  of  the  SCOUR  code.  In 
addition  to  providing  enhanced  insight  of  surface  burst 
dust  cloud  phenomena,  the  precursor  calculation  demonstrated 
the  feasibility  of  using  the  SCOUR  code  in  more  complicated 
flow  field  calculations. 

Chapter  IV  is  a presentation  of  SCOUR  results  for  the 
MIDDLE  GUST  (100  tons  of  TNT)  and  MIXED  COMPANY  (500  tons  of 
TNT)  high  explosive  detonations.  Near-surface  dust  cloud 
observations  and  concentration  measurements  (at  the  25  psi 
overpressure  station)  are  compared  to  SCOUR  predictions. 

The  earlier  SCOUR  results,  some  of  which  have  been  reported 
in  Reference  2,  were  very  encouraging  but  the  coding  errors 
(referred  to  above)  necessitated  a recalculation  of  the  HE 
dust  storms.  These  later  results  (SCOUR2-1975)  are  published 
in  Chapter  IV  to  provide  a convenient  comparative  study  with 
the  SCOUR  calculations  of  1973-74.  It  is  shown  that  the 
SCOUR2  predictions  differ  in  detail  from  the  earlier  results. 
For  large  diameter  pebbles  and  sand  grains,  the  aerodyncimic 
shear  lofting  is  greater  than  with  more  claylike  soils.  It 
was  also  shown  that  the  particle  impact  model  seriously  af- 
fects the  late  time  dust  concentrations  just  above  the  sur- 
face. 

The  fifth  and  final  chapter  is  a presentation  of  the 
S^  study  of  "reverse  percolation".  This  is  the  condition 
that  occurs  when  the  air  in  porous  soils  is  trapped  and 
pressurized  during  the  early  positive  phase  of  the  blast 
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wave  flow  field.  Initially,  the  air  flow  is  into  the  ground, 
but  when  the  overpressure  decays,  the  trapped  air  rapidly  per- 
colates up  through  the  soil  matrix.  Under  nuclear  blast  wave 
conditions  this  percolation  could  be  so  intense  as  to  blow  off 
surface  dust  and/or  "prepare"  the  surface  for  aerodynamic 
shear  lofting  of  soil  particles.  This  effect  has  been  ob- 

f 8 1 

served  in  experiments  conducted  in  an  inverted  shock  tube, ^ ^ 
and  was  analyzed  by  with  the  aid  of  a simple  model for 
flow-through  porous  media  in  Reference  2.  In  Chapter  V,  some 
examples  of  reverse  percolation  have  been  computed  with  the 
aid  of  a new  code  (FLUB)  that  was  developed  to  study  the 
dynamics  of  fluidized  beds.^^®'^^^  The  FLUB  calculations 
demonstrate  that  the  potential  exists  for  massive  blowoff  due 
to  reverse  percolation;  but  under  the  assumption  of  uniform 
fluidization,  the  velocities  are  too  small  for  particles  to 
be  actually  ejected  into  the  flow  field.  Rather,  the  effect 
would  be  to  make  particles  more  susceptible  to  aerodynamic 
shear  pickup  and  to  modify  the  impact  lofting  processes  cur- 
rently included  in  the  SCOUR  model. 


II.  SCOUR  CODE 
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The  SCOUR  code  has  been  described  in  the  other 
annual  reports.  In  this  chapter,  the  1973  version  of  SCOUR 
is  briefly  described  with  special  emphasis  on  those  sub- 
routines and  models  that  are  discussed  in  subsequent  chapters. 
An  appendix  is  provided  that  summarizes  1975  SC0UR2  nomencla- 
ture, input,  output,  and  special  diagnostic  routines. 

2.1  BASIC  DESCRIPTION 

SCOUR  was  originally  written  to  calculate  the  dust 
sweepup  contributions  to  the  nuclear  dust  cloud.  The  surface 
around  the  blast  center  is  divided  into  annular  rings.  Given 
an  axisymmetric  hydrodynamic  flow  field  input  (from  SHELL  or 
HULL  calculations) , SCOUR  calculates  the  amount  of  dust  mass 
lofted  from  each  ring,  the  time  at  which  it  is  ejected  from 
the  surface,  and  the  trajectories  of  representative  dust 
particles  as  they  move  in  the  near-surface  boundary  layer  or 
the  (free  stream,  hydrodynamic)  blast  wave  flow  field. 

The  air  equation  of  state,  spherical  particle  drag 
coefficient  routines  (including  Mach  number  dependence) , and 
the  routines  for  retrieving  the  hydrodynamic  flow  field  prop- 
erties for  the  trajectory  calculations  have  been  obtained 
from  the  AFWL's  DUSTY  code.  During  this  report  period,  a new 
subroutine  was  written  that  allows  tracer  particles  to  be  pre- 
loaded  at  specified  positions  and  velocities. 

A typical  SCOUR  calculation  proceeds  as  follows.  The 
preload  particles  are  input,  along  with  soil  particle  size 
distribution,  and  any  special  parameter  changes  from  the 
nominal  SCOUR  prescription.  The  radial  increment  for  the 
SCOUR  zones  are  scaled  by  the  cube  root  of  the  yield  of  the 
burst  that  was  the  source  for  the  input  hydrodynamic  flow 
field.  At  1 MT,  these  zones  are  50  meters  wide.  The  air 
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shock  is  located  in  the  hydro  flow  field  at  the  maximum  air 
density  position.  As  the  air  shock  traverses  a surface  zone, 
the  mass  of  dust  for  each  particle  size  to  be  lofted  via  aero- 
dynamic shfc^r  is  calculated  from  the  aerodynamic  mass  lofting 
equation  using  the  free  stream  air  flow  conditions.  This 
determines  the  mass,  M^,  to  be  represented  by  shear  lofted 
particles  from  that  surface  zone.  The  mass  lofting  rate  is 
then  time-integrated  for  each  subsequent  time  step  in  the 
calculation.  Whenever  a zone's  mass  integral  for  any  of  the 
particle  sizes  exceeds  M , or  some  specified  fraction  of  M , 

r tr 

a new  SCOUR  particle  is  created.  During  this  time  period  It 

was  decided  to  use  0.2  M as  the  mass  lofting  criteria  for 

P 

all  particles  launched  after  the  first  batch  of  shock  lofted 
particles  were  introduced  to  the  flow  field.  (This  has  the 
advantage  of  improving  the  resolution  of  the  SCOUR  code  by 
more  than  doubling  the  number  of  particles  to  be  tracked  in 
a SCOUR  calculation.)  During  the  positive  phase  it  is  pos- 
sible that  a finite  time  delay  exists  between  the  air  shock's 
arrival  and  significant  vertical  motion  of  the  activated  sur- 
face particles.  This  occurs  because  the  shear  first  moves 
particles  in  a nearly  horizontal  direction.  Should  the  delay 
time  for  the  particle  be  less  than  the  time  step  in  the  hydro 
flow  field  input,  the  representative  particle  is  ejected  into 
the  boundary  layer.  If  the  reverse  is  true,  particles  are 
stored  until  the  accumulated  time  exceeds  the  delay  time. 

All  particles  are  lofted  from  random  positions  in  the  zone. 

Particles  in  the  boundary  layer  are  subjected  to  tur- 
bulent diffusion  and  advection  by  the  local  air  flow.  Those 
particles  that  are  transported  out  of  the  boundary  layer  are 
moved  only  by  the  free  stream  wind  field.  As  the  winds  die 
down,  particles  may  impact  the  surface.  SCOUR' s impact  model 
enables  particles  with  sufficient  momentum  to  eject  new 
particles.  In  such  a manner,  t,||||^negative  phase  dust  storm 
will  degenerate  to  a saltation  type  of  dust  transport,  as 


described  in  Reference  12.  Should  a particle  impact  the 
surface  with  sufficient  momentum,  a new  particle  with  equal 
representative  mass  is  ejected  The  diameter  of  the  new 
particle  is  randomly  selected  from  the  in  situ  soil  distri- 
bution. The  mass  of  impacting  dust  is  subtracted  from  the 
aerodynamic  mass  lofting  integral  so  that  at  late  times 
there  is  very  little  mass  added  to  the  dust  cloud.  Finally, 
when  impacting  particles  do  not  have  sufficient  momentum, 
new  particles  are  not  ejected  and  the  dust  is  deposited. 

This  is  referred  to  as  a particle  "kill"  and  the  dust  mass 
is  subtracted  from  the  total  lofted  dust  mass  for  the  zone 
where  deposition  occurs. 

SC0UR2  has  been  initially  written  to  use  four  particle 
sizes  to  represent  the  in  situ  soil.  The  mass  distribution 
among  these  particles  is  prescribed  as  code  input.  (This  is 
not  mandatory,  however,  and  randomly  selected  partxcle  dia- 
meters between  certain  ranges  could  easily  be  .'.ncluded,  al- 
though the  price  of  each  calculation  would  increase  with  each 
additional  particle  size  class.) 

In  a typical  calculation,  approximately  5,000  particles 
are  tracked  in  SCOUR  at  a given  time.  Once  introduced  to  the 
flow  field,  particle  trajectories  are  calculated  during  hydro 
flow  field  time  steps.  Each  particle's  mass,  position,  and 
velocity  are  stored  as  it  is  moved  by  the  flow  field.  A 180 
second  surface  burst  SCOUR2  calculation  will  require  approxi- 
mately one  hour  of  computer  time  on  the  Univac  1108.  With  a 
complete  set  of  output,  a typical  price  per  calculation  is 
approximately  $300. 

The  on-line  output  of  SCOUR  is  geared  for  rapid  diac;  - 
nosis  of  the  dust  environment.  It  includes: 

• Particles'  mass,  position,  and  velocity  stored 
for  each  hydro  cycle  — available  for  printout. 

This  includes  the  tracking  of  pre-selected 


particles  (e.g.,  positive  phase  dust,  pre-loaded 
dust,  etc. ) . 

• Local  dust  concentrations  and  dust  dynamic  pres- 
sure for  combined  and  individual  particle  sizes 
in  a grid  superposed  on  the  dust  flow. 

• Total  mass  airborne. 

• Mass  above  given  surface  zone. 

• Mass  in  disc  between  two  altitudes. 

• Boundary  layer  details. 

• Local  erosion/deposition  depth. 

• Maximum  altitudes  of  particles. 

A second  off-line  output  format  is  also  available.  The 
code  user  specifies  the  overpressure  locations,  vertical  zones 
of  interest,  and  dust  cloud  parameters  of  interest.  The  new 
output  routine  averages  the  time  histories  of  the  specified 
variables  over  a spread  of  three  radial  zones  for  each  radial/ 
vertical  location.  The  spatial  averaging  eliminates  trouble- 
some "zeros"  in  the  digital  output.  Computer  plots  of  these 
time  histories  can  be  provided  per  user  requirements.  Cur- 
rently, SC0UR2  off-line  output  covers  the  following  variables: 

• Dust  concentration. 


• Dust  flux,  horizontal. 

• Dust  flux  impulse. 

• Dust  flux,  vertical. 

• Dust  flux  impulse,  vertical. 

• Mass  averaged  dust  velocity,  horizontal. 

• Mass  averaged  dust  velocity,  vertical. 

• Dust  dynamic  pressure. 

• Dust  dynamic  pressure  impulse. 
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The  input-output  format  for  SC0UR2  is  fully  described  in  the 
appendix  to  this  report. 

2.2  NON- INTERACTIVE  FORMULATION 

The  basic  SCOUR  model  was  formulated  during  the  first 
year  of  the  dust  sweepup  program.  A number  of  basic  assump- 
tions were  made  at  the  outset  that  define  the  scope  of  the 
modeling  problem.  These  were  based  on  measurements  of  dust 
clouds  and  earlier  calculations  of  the  dust  cloud which 
indicated  that  only  in  the  near-surface  region  may  particulate 
concentrations  become  dense  enough  to  measurably  affect  the 

* ri61 

local  flow  field.  Typical  dust  loading  ratios  indicate 

that  dust  concentrations  in  the  stem  and  mushroom  portion  of 

— 6 — 1 0 

the  cloud  are  between  10  g/cc  and  10  g/cc.  It  is  entirely 
reasonable,  therefore,  to  approach  the  dust  lofting  problem 
with  the  following  non- interactive  formulation; 

• Particles  are  moved  by  the  flow  field  but  do  not 
affect  it. 

• Characteristic  particles,  representing  a large 
number  of  real  particles,  are  introduced  to  the 
flow  field  in  lieu  of  following  each  and  every 
particle  or  treating  the  raised  dust  statistically. 

• Trajectories  of  characteristic  particles  are  inde- 
pendent of  other  particles. 

• The  boundary  layer  at  the  surface  may  be  computed 
under  the  assumption  that  the  input  hydrodynamic 
flow  field  provides  the  free  stream  information. 

The  non-interactive  model  is  an  important  simplifica- 
tion as  it  precludes  the  expenditure  of  time,  funds,  and  ef- 
fort on  reformulating  the  air  flow  hydro  codes.  Hence  for  a 

* 

This  was  confirmed  in  the  original  (and  subsequent)  SCOUR  cal- 
culations. 
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given  weapon  yield  and  HOB,  the  most  appropriate  calculation 
in  the  tape  library  of  SHELL  or  HULL  hydro  calculations  can 
be  used  to  compute  the  dust  cloud. 

2.3  SC0UR2  BOUNDARY  LAYER  MODEL 

The  boundary  layer  model  in  SC0UR2  is  almost  identical 
to  that  described  in  References  1 and  2.  It  was  derived  on 
the  basis  of  a set  <^f  simplifying  assumptions,  consistent  with 
the  non-interactive  formulation,  that  make  the  problem  tract- 
able within  the  scope  of  the  overall  study.  One  of  these  is 
the  assumption  that  the  roughness  effect  of  the  dust  particles 
moving  just  above  the  surface  is  similar  to  that  produced  by 
very  rough  walls,  i.e.,  the  "perfectly  rough"  regime  of  turbu- 
lent boundary  layers.  The  shear  in  the  near-surface  "perfectly 
rough"  boundary  layer  is  dependent  only  on  free  stream  velocity 

and  a suitably  selected  "effective  roughness".  On  the  basis  of 

fl71 

high  velocity  wind  tunnel  tests,  it  was  determined  that  the 

average  height  of  the  near-surface  particles  provides  an  ade- 
quate measure  of  effective  roughness. 

The  SCOUR  formulation  has 


u* 

U 


0.707 


'2.87  + 0.7  Jin 


y' 


(2.1) 


where 


= shear  velocity  in  the  boundary  layer 
= /p  , T = shear  at  wall 

V 0 0 


P = 


0 0 
air  density 


U = free  stream  velocity 


distance  behind  air  shock  for  the  positive  phase, 
and  y*U  dt  in  the  negative  phase  (t  is  time) 
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kg  = roughness  height 

The  maximum  ratio  of  u*/U  is  0.248,  but  in  our  calculations 

^ e 

0.15  > ^ > 0.03  is  a typical  range  of  values. 

From  the  "perfectly  rough"  shear  law,  it  was  also  pos- 
sible to  derive  integral  expressions  for  the  blast  wave  boundary 
layer  that  are  evaluated  under  the  assumptions  that  the  local 
air  flow  is  quasi-steady  (i.e.,  keyed  to  shock  motion),  incom- 
pressible, and  the  horizontal  velocity  in  the  boundary  layer 
is  governed  by 

/ \l/n 

S-  = (^)  ,y<«  (2.2) 


where 

u = radial  component  of  air  velocity 

Y = altitude  above  original  ground  surface 

6 = boundary  layer  thickness 

n = boundary  layer  profile  exponent,  8 > n ^ 2; 

(n  = 2 in  SCOUR,  on  the  basis  of  wind  tunnel  ex- 

[171 

periments  ) . 

The  turbulent  boundary  layer  development  behind  the 
blast  wave  from  a 100  ton  HE  detonation  (MIDDLE  GUST)  was 
measured  at  the  30  and  70  psi  overpressure  locations. 

These  results  for  the  blast  wave  boundary  layer  and  the  pre- 
dictions of  turbulent  boundary  layer  growth  behind  the  air 
shock  are  provided  by  more  sophisticated  boundary  layer  models 
and  compared  to  the  SCOUR  boundary  layer  predictions  in  Fig- 
ure 2.1.  The  SCOUR  results  at  early  times  are  independent  of 
particle  roughness  effects  because  the  topographical  rough- 
ness dominates  the  SCOUR  boundary  layer  development.  It  is 
interesting  to  note  that  the  SCOUR  results  are  in  good  agree- 
ment with  the  other  theories,  except  at  later  times  where  the 
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Figure  2.1a.  Comparison  of  SCOUR  boundary  layer  thickness  to 
other  predictions  and  the  experimental  data  at 
the  nominal  70  psi  overpressure  location  (280  ft 
station)  for  MIDDLE  GUST.  [18] 
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Figure  2.1b. 


Comparison  of  SCOUR  boundary  layer  thickness  to 
other  predictions  and  the  experimental  data  at 
the  nominal  30  psi  overpressure  location  (260  ft 
station)  for  MIDDLE  GUST.  [18] 


rough  wall  assumption  leads  to  significantly  greater  boundary 
layer  growth.  The  major  reasons  for  the  slight  underpredic- 
tion at  early  times  are  that  the  shear  law  (Eq.  (2.1))  has  a 
maximum  shear  limit  that  reduces  the  initial  growth  of  the 
boundary  layer  and  that  code  resolution  at  early  times  is  in- 
sufficient to  accurately  compute  initial  growth  even  when  the 
shear  is  below  the  upper  limit. 

In  any  SC0UR2  calculation,  there  is  a regime  of  non- 

uniform  behavior  and  rapid  boundary  layer  growth  associated 

with  the  transition  from  the  positive  to  negative  phase. 

Generally,  the  velocity  and  density  gradients  are  relatively 

small  during  flow  reversal.  An  exact  formulation  of  the 

problem  would  have  to  treat  the  interaction  between  the 

boundary  layer  and  the  free  stream.  To  avoid  this  difficult 

(but  not  crucial)  portion  of  the  boundary  layer  calculation, 

a maximum  boundary  layer  height,  DELMAX,  is  prescribed  for 

each  SCOUR  calculation.  DELMAX  was  originally  based  on  the 

maximum  size  of  the  zones  in  the  hydrodynamic  flow  field  to 

be  used  in  the  SCOUR  calculation.  It  was  later  determined 

that  SCOUR  predictions  were  insensitive  to  realistic  changes 
r 2 1 

in  DELMAX.  Currently  the  prescription  for  DELMAX  is 

DELMAX  = 175  (W)^"^^  (2.3) 

where  W is  the  yield  in  MT. 

It  should  be  noted  that  a minor  programming  error  in 
the  subroutine  for  transition  to  the  negative  phase  boundary 
layer  calculation  was  noted  during  this  report  period.  The 
error  had  the  effect  of  delaying  the  negative  phase  calcula- 
tion by  one  hydro  time  step.  It  was  corrected  at  the  end  of 
the  report  period. 
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2.4 


AERODYNAMIC  LOFTING 


The  aerodynamic  mass  lofting  prescription  in  SC0UR2  is 
based  on  a simple  premise. 

The  force  required  to  loft  particles  out  of  the 
fluidized  surface  layer  is  supplied  by  some  fraction, 

layer. 


of  the  net  shear  in  the  near-surface  boundary 


For  a given  particle  size. 


M V = b„^ 
p 0 ND 


Pa'^*  - ^ 


th  ‘^i 


(2.4) 


where 


ND 


= shear  transfer  coefficient 


= threshold  shear  (taken  as  1 dyne/cm^ ) 

dA^  = surface  area  of  zone 

M = lofting  rate  (at  dA) 

P 

V = initial  ejection  speed  of  particulates 
0 

Conceptually,  b^j^,  and  are  all  dependent  on; 

particle  size,  soil  type,  local  wind  conditions,  soil  surface 
cohesion,  vegetation  cover,  etc.  In  SC0UR2,  the  nominal  pre- 
scription is  that  there  is  an  overall  10  percent  conversion  of 

near-surface  aerodynamic  shear  to  the  vertical  force  trans- 

★ 

mitted  to  lofted  particles.  The  lofting  efficiency  for  each 
particle  size  class  is  directly  proportional  to  its  mass 
fraction,  i.e. , 


^ND  ^D^N 


(2.5) 


* (221 

According  to  Bagnold,  '■  ^b.,  could  be  as  high  as  0.3  for  dry, 

hard  soils  that  are  characterized  by  both  large  and  small 
surface  particles. 
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where 


= mass  fraction  of  soil  represented  by  particle 
with  diameter,  D 

= overall  shear  conversion  coefficient  (nominally 

0.1) 

The  remainder  of  the  shear  is  assumed  to  be  dissipated 
in  the  near-surface  tangential  motion  of  the  particulates  (the 
so-called  "creep"  layer) . Hence,  the  SCOUR  code  tracks  only 
those  particles  which  are  ejected  out  of  the  near-surface 
"creep"  layer. 

(21  231 

Based  on  recent  observations  ’ of  particle  entrain- 
ment at  high  velocities,  the  particle  ejection  speed,  V , is 

0 

randomly  selected  between  0.5  u^  and  u*,  with  initial  direction 
chosen  so  that  the  slowest  particles  are  ejected  at  an  angle 
of  10°  to  the  horizontal  while  the  fastest  are  ejected  at  90°. 
The  horizontal  component  of  the  particle  velocity  was  pro- 
grammed to  be  in  the  direction  of  the  positive  phase  flow, 
even  in  the  negative  phase.  This  error  was  detected  at  the 
end  of  this  report  period  and  the  present  version  of  SCOUR2 
contains  the  prescription  that  ejects  particles  with  horizontal 
velocities  in  the  same  direction  as  the  instantaneous  local 
air  flow  above  a surface  zone. 

2.5  IMPACT  LOFTING  AND  DEPOSITION 

When  the  wind  speeds  decay,  the  bigger  particles  will 
fall  back  to  the  ground  (particularly  those  above  1 mm  in  dia- 
meter) . The  SCOUR  impact  model  treats  this  situation  by  re- 
ejecting a particle,  whose  diameter  is  randomly  selected  from 
the  in  situ  particle  size  distribution,  and  which  represents 
the  same  mass  as  the  impacting  particle.  It  is  ejected  verti- 
cally with  0.1  of  the  horizontal  momentum  and  all  of  the  verti- 
cal momentum  of  the  impacting  particle.  The  mass  ejected  is 


sur- 


L 


subtracted  from  the  mass  lofting  integral,  for  the  i^ 
face  zone. 


(mp) 


AA. 

1 


/Apdt, 


(mass  lofting  integral) 


(2.6) 


which  computes  the  amount  of  mass  being  lofted  due  to  shear 
for  a given  particle  size  (M  from  Eq.  (2.5)).  Hence,  as 
particles  fall  out,  the  mass  lofting  integrals  at  a given 
surface  zone  may  no  longer  increase,  resulting  in  the  estab- 
lishment of  the  equilibrium  saltation  mode  of  particulate 
motion  at  the  surface  (wherein  no  net  erosion  or  deposition 
is  occurring) . 

Experiments  in  both  natural  and  laboratory  conditions 
have  revealed  that  there  are  threshold  levels  of  shear  in  the 
boundary  layer,  below  which  saltation  ceases. In  the 
shear  threshold  regime,  impacting  particles  begin  to  deposit 
(rather  than  cause  the  ejection  of  other  particles) . As  one 
might  expect,  this  process  is  dependent  on  soil  type  and 
particle  size.  Figure  2.2  shows  the  variation  of  the  thres- 
hold shear  velocity,  with  particle  size  for  uniform 

sandy  soils. 

The  particle  deposition  model  in  SCOUR  is  based  on  the 
principle  that  particle  motion  at  the  surface  is  directly 
linked  to  the  existence  of  the  threshold  shear.  That  is, 
particle  impacts  will  no  longer  be  self-sustaining  if  the 
ejection  velocity  for  a SCOUR  particle  is  less  than  the 
threshold  ejection  velocity  in  Figure  2.2).  Thus,  for 

an  impact  where  the  particle  ejection  velocity  is  less  than 
u*th'  SCOUR  "kills"  the  impact  particle  and  the  deposition  is 
recorded  in  the  surface  zone  erosion  integral.  When  subse- 
quent new  particles  are  created  by  the  action  of  aerodynamic 
shear,  they  replace  the  killed  particles  in  the  SCOUR  data 
storage  subroutine. 
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Values  of  u.. 


Figure  2.2.  Plot  of  the  fluid  threshold  and  impact  threshold 
shear  velocities  that  are  necessary  to  induce 
saltation  (from  Reference  22).  For  aerodynamic 
lofting,  the  fluid  threshold  is  appropriate 
= 1.0  dynes/cm  used  in  SCOUR).  Impact 
lofting  is  governed  by  the  impact  threshold. 


2.6 


PARTICLE  ENTRAINMENT 


The  movement  of  the  enormous  quantities  of  dust  raised 
by  the  blast  wave  winds  is  treated  in  SCOUR  by  monitoring  the 
motion  of  numerous  5,000)  individual  particles,  each  repre- 
senting a much  larger  number  of  particles  in  a certain  diameter 
class.  The  trajectories  of  the  representative  particles  are 
computed  by  integrating  the  equations  of  motion  in  time  steps 
corresponding  to  the  time  sequence  of  data  "dumps"  in  the 
hydro  flow  field  being  used  in  a particular  SCOUR  calculation. 
The  equations  are 


m 


P 


m 


P 


-"p9  * Cd(I  Pa)  (''p-''a)  (^p-''a) 

s(l  Pa)  (''p-'^a)  (^p-“a) 


(2.7) 


(2.8) 


where 

Vp  = speed  of  particle  = 

= local  wind  speed  = ^u^  + v^ 

• • 

Y , X = components  of  particle  velocity  in  y,  x 

directions 

V , u = component  of  wind  velocity  in  y,  x directions 

Ci  a 

The  drag  coefficient,  C^,  is  a function  of  Reynolds  number  (Re) 
and  Mach  number  (both  based  on  the  relative  speed  between  air 
and  particle)  . A subroutine  for  C^^  and  the  air  equation  of 
state  was  extracted  from  AFWL's  DUSTY  code. 

2.6.1  Diffusion  Model 

It  is  evident  from  Eqs.  (2.7)  and  (2.8)  that,  except  for 
the  initial  conditions  given  to  each  particle,  all  dust  motion 
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is  governed  by  the  action  of  aerodynamic  drag  and  gravity 
forces  exerted  on  the  individual  dust  particles.  The  hydro 
flow  fields  used  as  input  to  the  SCOUR  calculations  do  not 
include  boundary  layer  effects  or  turbulent  diffusion  in  the 
free  stream.  Hence,  particles  injected  into  an  unmodified 
flow  field  can  only  be  entrained  and  suspended  if  the  vertical 
components  of  the  air  velocity  are  great  enough. 

SCOUR  takes  into  account  the  diffusion  of  particles  in 
the  boundary  layer  by  prescribing  diffusional  (or  transport) 
velocities  that  are  superposed  on  the  horizontal  and  vertical 
velocity  boundary  layer  profiles.  Hence,  particles  will  be 
suspended  and/or  transported  out  of  the  boundary  layer  depend- 
ing on  the  particles'  susceptibility  to  the  prescribed  trans- 
port velocity. 

The  SCOUR  model  for  transport  velocity  is  based  on  the 
analogy  to  diffusion  in  atmospheric  boundary  layers.  Assuming 
that  dust  concentrations  will  change  over  the  same  spatial 
scale  as  the  dominant  mixing  eddies,  it  can  be  shown  that  the 
transport  velocity  in  the  boundary  layer,  V^,  is  approximated 
by 


V-  = (0.45)  a^u 
^ 0 


(2.9) 


where 

a = standard  deviation  of  the  horizontal  wind  angle 
e 

fluctuation 

u = local  horizontal  wind  velocity 

This  suggests  that  a simple  prescription  for  would  be  of 
the  form 


where 


-k'  = 0.45  . (2.11) 

This  is  the  diffusion  prescription  in  the  SCOUR2  model. 

The  observed  values  of  a yields  values  of  k'  in  the 

[2]  ® 

range 

0.004  < k"  < 0.12  , (2.12) 

where  the  lower  end  of  the  range  is  associated  with  very  stable 
stratified  boundary  layers  (i.e.,  cold  air  near  the  ground, 
warmer  air  above) . At  neutral  stability  where  the  thermal 
gradient  matches  the  temperature  loss  associated  with  the 
adiabatic  rise  of  a gas  particle  through  the  atmosphere, 
k'  » 0.06.  The  high  end  of  the  spectrum  is  associated  with 
unstable  atmospheres. 

The  nominal  value  for  k'  is  0.08.  This  value  gives  good 
agreement  for  the  dust  cloud  concentration  measurements  ob- 
tained from  the  MIDDLE  GUST  and  MIXED  COMPANY  HE  detonations. 

It  could  be  anticipated  that  for  the  case  of  a preheated  sur- 
face layer,  the  dust  transport  would  be  intensified  due  to 
boundary  layer  instabilities.  For  this  reason,  k""  is  0.20 
for  precursor  cloud  calculations.  The  advection  of  the  dust 
behind  the  shock  wave  in  these  cases  is  dominated  by  the  free 
stream  precursor  flow  field,  and  the  exact  value  of  k'  does 
not  play  a significant  role  in  dust  cloud  development. 

To  account  for  the  effects  of  near-surface  diffusion 
phenomena  in  the  SCOUR2  model,  numerical  turbulence  was  super- 
posed on  the  boundary  layer  velocity  prescription  to  simulate 

(21 

a randomized  mixing  process  due  to  local  disturbances. 
Specifically, 

= k'u  + 0.1  u^,  y £ 6 (2.13) 

where  the  sign  in  front  of  0.1  u^  is  randomly  chosen. 
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2.6.2  Entrainment  at  the  Shock 

The  SCOUR  boundary  layer  model  does  not  compute  the 
boundary  layer  parameters  above  a surface  zone  until  the  air 
shock  has  completely  traversed  the  zone.  As  a result,  parti- 
cles which  are  being  swept  along  by  the  shock  and  are  in  the 
same  radial  zone  as  the  shock  will  not  be  subjected  to  either 
of  the  entrainment  models  discussed  above.  To  compensate  for 

the  enhanced  entrainment  and  diffusion  that  has  been  observed 

[21] 

in  shock  tube  experiments,  the  velocity  field  in  the  shock 

radial  zone  is  prescribed  as  follows: 


u = U 


i-k) 


where 


V = 0. 08  U + U ( r; i 

e e \ X„  - X . / 

' S IS  / 


= particle  altitude 


(2.14) 


(2.15) 


^is'  ^ 


= radial  location  of  shock 

= radial  location  and  boundary  layer  height  at 
nearest  zone  completely  traversed  by  shock 


III.  NUCLEAR  DUST  CLOUDS 

Two  sets  of  dust  cloud  calculations  were  made  during 
this  report  period.  The  first  group  (I)  consisted  of  cal- 
culations for  32.5  KT  and  5 MT  surface  bursts  and  a 36.6  KT 
detonation  at  350  ft  HOB.  These  were  primarily  intended  to 
provide  "state  of  the  art"  dust  cloud  information  and  order 

of  magnitude  comparative  predictions  for  input  to  other 

r 5-71 

defense  contractors,  prior  to  the  setting  of  nominal 

values  for  SCOUR  parameters.  These  results  are  summarized 
in  Sections  3.1  and  3.2. 

The  second  set  of  nuclear  dust  cloud  calculations  (II) 

were  made  for  surface  burst  yields  of  0.1  MT , 1.0  MT,  and 

20  MT.  SCOUR  parameters  were  set  at  nominal  values  so  that 

the  effect  of  scale  on  the  SCOUR  dust  cloud  predictions  could 

be  studied.  The  particle  size  distribution  of  the  soil  was 

the  same  as  that  used  in  most  of  the  previous  SCOUR  calcula- 
f 1 21 

tions;  ' i.e.,  representative  particles  of  0.8  cm,  0.3  cm, 

and  0.07  cm  in  diameter  were  assigned  20  percent  each  of  the 
total  soil  mass.  The  remaining  40  percent  was  assigned  to 
representative  particles  with  a diameter  of  0.005  cm.  These 
results  are  summarized  in  Section  3.3. 

The  SCOUR  code  results  presented  in  this  chapter,  as 
has  been  previously  explained,  have  been  superseded  by  later 
calculations  published  in  References  3 and  4.  Nevertheless, 
they  provide  a qualitative  indication  of  the  SCOUR  code's  pre- 
dictive capacity. 

3.1  SURFACE  BURST  DUST  CLOUD  PARAMETERS  - I 

A number  of  SCOUR  calculations  were  made  in  the  first 
half  of  the  report  period,  at  yields  between  32.5  KT  and 
5.0  MT.  For  the  purposes  of  this  report,  the  results  for  the 
32.5  KT  and  5.0  MT  surface  bursts  are  briefly  reviewed  to 
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illustrate  the  sensitivity  of  the  code  to  changes  in  the 

f 1 21 

diffusion  model  in  the  old  ' and  newer  version  of  SCOUR 
(fully  described  in  Reference  2).  The  32.5  KT  calculation 
was  made  with  a scaled  version  of  AFWL/SHELL  flow  field 
No.  1002.055  (100  KT  surface  burst,  0-300  secs).  The  5.0  MT 
calculation  utilized  AFWL  problem  No.  143.000  (1.4  MT  sur- 
face burst,  0-600  secs) . Both  used  the  nominal  in  situ 
particle  distribution  (20  percent  mass  to  0.8  cm,  0.3  cm, 

0.07  cm  diameter,  40  percent  mass  to  0.005  cm  diameter). 

Typical  concentrations  in  the  dust  sweepup  cloud  for 
a 5 MT  surface  burst  are  presented  in  Figure  3.1.  During 

-3 

the  positive  phase,  these  range  from  2 x 10  g/cc  to 

_ 3 

0.2  X 10  g/cc,  over  an  overpressure  range  of  500  psi  to 
2 psi.  Negative  phase  concentrations  are  predicted  to  be  ap- 
proximately 0.25  X 10  ^ g/cc  with  the  old  SCOUR  diffusion 
model, and  0.05  x 10~^  g/cc  in  the  (more  diffusive)  new 
SCOUR  calculations.  (This  is  the  present  diffusion  model  in 
SCOUR.)  Predicted  values  of  the  near-surface  cloud  height 
as  a function  of  radial  distance  from  ground  zero  (at  t = 

180  secs)  are  plotted  in  Figure  3.2  for  a (typical)  surface 
burst. 

The  mass  loading  of  the  cloud  for  two  new  SCOUR  cal- 
culations for  a 32.5  KT  surface  burst  are  presented  in  Fig- 
ure 3.3.  These  curves  represent  mass  (MT)  per  MT  yield  that 
is  lofted  above  the  indicated  altitudes.  In  one  case,  the 
diffusion  constant,  (in  Eq.  (2.10)),  was  maintained  at 
0.225.  These  results  are  compared  to  that  of  a second  calcu- 
lation wherein  the  diffusion  velocities  were  reduced  to  one- 
tenth  the  previous  values  (k^  = 0.0225). 

The  order  of  magnitude  change  in  k'',  the  diffusion 
parameter,  has  the  greatest  effect  on  sweepup  contribution  to 
the  high  altitude  dust  cloud.  These  results  are  similar  to 
that  obtained  in  the  1 MT  calculations  presented  in  Figure  6.7 
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Figure  3.2.  Variations  of  cloud  height  with  range  for  a 5 MT 

surface  burst.  Points  represent  altitude  of  cells 
above  individual  surface  zones  populated  with  dust 
above  lO”  g/cc.  Note  that  these  results,  while 
qualitatively  consistent  with  subsequent  SC0UR2 
calculations,  are  superseded  by  results  published 
in  References  3 and  4. 
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Figure  3.3. 


Dust  loading  due  to  sweepup.  Curves  show  mass 
lofted  above  ground  level  (0  m)  and  300  m for  a 
32.5  KT  (surface  burst)  using  two  different 
values  of  k'  (kj  = 0.255).  Note  that  these  re- 
sults, while  qualitatively  consistent  with  subse 
quent  SCOUR2  calculations,  are  superseded  by  re- 
sults published  in  References  3 and  4. 
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of  Reference  2.  In  the  present  calculations,  the  cloud  rise 
velocity  is  57  m/sec  (compared  to  115  m/sec  for  the  1 MT 
case  in  Reference  2) . It  was  concluded,  therefore,  that  the 
50  percent  variation  in  total  cloud  loading,  and  order  of 
magnitude  variation  in  high  cloud  loading,  represents  the 
general  sensitivity  of  the  SCOUR  code  to  the  possible  range 
of  values  for  the  diffusion  parameter. 

Figure  3.4  is  a plot  of  the  predicted  depth  of  soil 
eroded  from  the  surface  at  t = 18.3  secs  and  t = 168  secs 
for  a 5 MT  surface  burst.  Peak  erosion  (1.35  cm)  due  to 
aerodynamic  shear  occurs  at  about  2 km,  the  70  psi  air  shock 
overpressure  location.  About  40  percent  of  the  erosion  oc- 
curs during  the  positive  phase.  The  appearance  of  a maximiim 
in  erosion  at  R ~ 2 km  (rather  than  a monotonic  increas.e  in 
erosion  as  ground  zero  is  approached)  is  due  to  smaller  air 
dynamic  pressure  impulse  in  the  positive  phase  for  R < 2 km, 
and  relatively  smaller  horizontal  velocities  at  close-in 
locations  during  the  negative  phase  (corresponding  to  the 
stem  region  containing  large  updraft  velocities) . 

3.2  LOW  HEIGHT  OF  BURST  DUST  CLOUD  - 36.6  KT  AT  350  FT  HOB 

The  dust  clouds  raised  by  the  air  blast  waves  associated 
with  the  various  tower  shots  at  the  Nevada  Test  Site  (NTS)  have 
been  recorded  on  film.  These  motion  picture  films  provide  a 
reservoir  of  data  to  make  qualitative  comparisons  between  the 
calculated  and  observed  dust  clouds  created  by  non-cratering 
nuclear  explosions.  During  this  report  period  an  important 
review  of  the  NTS  dust  cloud  data  was  initiated. The 
focus  of  this  study  was  to  develop  a phenomenological  under- 
standing and  model  for  the  development  of  precursors  to  the 
air  blast  waves.  As  described  in  Reference  25,  the  precursors 
are  distortions  in  the  blastwave  structure  caused  by  the 
radiant  preheating  of  an  air  layer  just  above  the  ground  sur- 
face. 


30 


! j 

Systems,  Science  and  Software  obtained  (and  converted  j ! 

for  use  with  SCOUR)  the  hydro  output  of  an  AFWL  SHELL  calcula-  ! 

tion  of  a precursor  flowfield  for  a 36.6  KT  yield  at  350  ft  j 

HOB  (AFWL  No.  362.017,  0.09  sec-5.0  secs).  In  this  AFWL 

I problem,  the  air  near  the  ground  was  preheated  prior  to  the 

inception  of  the  calculation.  The  preheat  layer  was  roughly 

I 5 meters  high  and  extended  out  to  a range  of  420  meters.  The 

presence  of  a preheated  layer  causes  the  shockwave  structure 
! ' 
I to  "toe  out"  just  above  the  ground  surface,  leading  to  higher 

I 

f dynamic  pressures  and  lower  overpressures  when  compared  to 

f;  ideal  surface  flow  fields. 

Preheating  was  achieved  by  designating  a fraction  of 
the  incident  theirmal  energy  on  the  ground  (from  the  fireball) 
to  be  transferred  to  the  air  layer  above  the  ground.  The 
flux  of  incident  energy  at  a slant  range,  r . , was  assumed 
to  be  equal  to  {eh/AiJr^^)  where  e is  the  fireball  flux  and 
h is  the  burst  height.  AFWL  also  assumed  that  the  fraction 
of  incident  energy  transferred  to  the  air  layer  was  equal  to 
the  ratio  of  h to  r.  To  simulate  the  popcorning  threshold 
' of  dust,  the  air  layer  was  not  heated  until  5.5  calories  per 

square  centimeter  (or  some  multiple  thereof)  were  accumulated. 

The  precursor  flow  fields  resulting  from  these  imposed 
heating  laws,  such  as  those  presented  in  Reference  27,  were 
. qualitatively  similar  to  observed  flow  fields  but  were  dif- 

ferent in  detail.  For  the  purposes  of  the  SCOUR  calculation, 
however,  we  were  most  interested  in  demonstrating  that  SCOUR 

l| 

could  be  used  with  multiple  shockwave  flow  fields  and  that  the 
resulting  dust  cloud  structure  was  in  qualitative  agreement 

I 

I with  observed  precursor  dust  clouds. 

Figure  3.5  is  a plot  of  the  air  velocity  vector  in  the  , 

I 

hydro  flow  field  for  two  times.  At  0.38  seconds,  the  precursor 


region  is  shown  to  consist  of  a distorted  velocity  field,  with 
high  updraft  velocities  just  behind  the  leading  shock.  These 
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thermal  perturbations  leave  a wake  of  local  eddy  disturbances 
in  the  negative  phase  flow  field,  but  the  leading  shock 
eventually  overtakes  the  precursor  and  a more  conventional 
shock  structure  develops  as  demonstrated  by  the  air  flow  at 
1.9  seconds  shown  in  Figure  3.5b. 

With  this  hydro  field  as  input,  a SCOUR  calculation 
of  the  dust  cloud  was  made.  The  early  positive  phase  results 
are  compared  to  a NTS  shot  of  similar  yield  and  height  of 
burst  (HARRY)  in  Figure  3.6.  These  results  indicate  that  the 
predicted  height  of  the  precursor  cloud,  which  moves  out  as 
hemi-torous  behind  the  shock,  is  quantitatively  and  qualita- 
tively comparable  to  the  NTS  observations.  An  additional  cal- 
culation using  a non-precursor  flow  field  (No.  362.018)  at  the 
same  yield,  gave  cloud  height  predictions  one-half  those  shown 
in  Figure  3.6.  This  indicates  that  the  precursor  effect  is 
important  to  include  in  dust  cloud  calculations  when  theory 
predicts  a precursor  will  be  evident. 

An  anomalous  result  in  this  preliminary  precursor  cal- 
culation, due  to  a coding  error  then  present  in  SCOUR  boundary 
layer  diffusion  models,  was  most  evident  in  the  late  time  dust 
cloud  predictions.  The  error  in  the  subroutine  lead  to  ex- 
tremely large  diffusion  velocities  as  particles  approached  the 
outer  edge  of  the  boundary  layer.  In  effect,  particles  would 
slowly  diffuse  through  the  boundary  layer  and  then,  inadvert- 
ently, be  ejected  out  of  the  boundary  layer.  This  effect  was 
not  immediately  recognized,  and  the  results  were  reported 
in  Reference  28  as  follows: 

"The  agreement  between  predicted  cloud  dimensions  and 
the  NTS  films  extend  to  the  early  development  of  the 
negative  phase  cloud.  Results  of  the  SCOUR  precursor 
calculations  are  compared  to  traces  of  the  TURK  event 
(43  KT  at  300  ft)  dust  cloud  taken  from  Reference  28. 

In  Figure  10  the  representative  particles  (dp  = 

0.005  cm)  entrained  in  the  flow  field  are  plotted  on 
the  same  axes  as  traces  of  the  cloud.  It  is  clearly 
evident  that  the  bulk  of  the  predicted  distribution 
of  material  lies  within  the  dust  cloud  traces. 


Altitude,  meters  Altitude,  meters 


Figure  3.6.  Comparison  of  predicted  precursor  cloud  heights 
(at  t = 0.23  and  0.7  sec)  to  that  obtained  from 
motion  pictures  of  the  HARRY  event. [24]  Note 
that  these  results,  while  qualitatively  consistent 
with  subsequent  SCOUR2  calculations,  are  super- 
seded by  results  published  in  References  3 and  4. 


This  result  lends  further  confidence  to  the  general 
boundary  layer-diffusion  model  in  SCOUR.  If  there 
were  no  diffusion,  the  near-surface  dust  cloud  would 
be  extremely  localized,  even  with  a precursor.  This 
is  demonstrated  in  Figure  11  where  the  minimal  cloud 
calculated  with  zero  diffusion  is  compared  to  the 
nominal  calculation  and  the  TURK  cloud.  It  is  evi- 
dent that  diffusion  is  required  in  the  model  to 
accurately  predict  cloud  dimensions." 

Once  the  diffusion  subroutine  error  was  recognized,  it 

became  obvious  that  the  late  time  cloud  profiles  in  Figures 

3.7  and  3.8  would  be  markedly  higher  than  those  predicted  by 
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a correct  version  of  the  code.  Subsequent  calculations  ' 
indicate  that  the  SCOUR  predictions,  without  the  free  stream 
diffusion  of  the  dust  particles,  actually  fall  in  between  the 
two  cloud  outlines  in  Figure  3.7.  Moreover,  the  puffs  of 
dust  at  very  high  altitudes  in  Figure  3.7  were  solely  the 
result  of  the  subroutine  error. 

These  problems  notwithstanding,  this  first  calculation 
with  a precursor  flow  field  did  provide  a demonstration  that 
SCOUR  could  be  utilized  (on  a first-pass  basis)  to  calculate 
the  precursor  dust  storm.  The  code  was  later  corrected  and 
modified  (in  this  report  period)  to  include  preload  particles 
that  represented  popcorn  dust.  No  further  precursor  calcula- 
tions were  made  in  this  report  period,  but  subsequent  SC0UR2 

r 3 4 ] 

calculations  ' have  provided  much  valuable  information  on 
precursor  dust  cloud  phenomenology. 

3.3  NUCLEAR  SURFACE  BURSTS  - (II) 

The  following  calculations  were  made  subsequent  to 
those  presented  in  Sections  3.1  and  3.2.  The  anomalous  dif- 
fusion at  the  outer  edge  of  the  boundary  layer  was  corrected, 
but  as  previously  mentioned,  some  minor  coding  errors  still 
were  present  in  SCOUR.  Nevertheless,  these  calculations  pro- 
vided very  close  estimates  of  the  corrected  SCOUR  results 
that  are  described  in  References  3 and  4.  The  diffusion 
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Figure  3.7.  Comparison  of  SCOUR  code  predictions  of  positions 
of  representative  dust  particles  to  traces  of 
cloud  produced  during  the  TURK [29]  event.  These 
results  are  superseded  by  those  published  in 
References  3 and  4.  They  are  qualitatively  in- 
consistent with  the  later  results  due  to  a coding 
error  that  magnified  the  diffusion  at  outer  edge 
of  the  boundary  layer. 
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Figure  3.8.  SCOUR  predictions  of  the  near-surface  dust  cloud 
height  with  and  without  diffusion,  at  5 seconds 
compared  to  the  TURK  cloud  at  6.5  seconds.  These 
results  are  superseded  by  those  published  in 
References  3 and  4.  They  are  qualitatively  in- 
consistent with  the  later  results  due  to  a coding 
error  that  magnif  d the  diffusion  at  outer  edge 
of  the  boundary  layer. 


parameter,  k' , was  0.08  for  both  the  positive  and  negative 
phases  of  these  calculations,  and  the  nominal  soil  particle 
size  distribution  was  assumed. 

3.3.1  Mass  Lofting  Capacity 

The  amount  of  dust  lofted  by  the  air  blast  wave  flow 
field  in  three  different  SCOUR2  calculations  is  plotted  in 
Figure  3.9.  Dust  mass  has  been  non-dimensionalized  by  the 
yield,  and  time  is  scaled  by  the  cube  root  of  the  yield  (in 
MT) . Figure  3.9  also  indicates  the  total  dust  mass  deposited 
and  the  dust  lofted  into  the  high  cloud  (Y  > 3 km) . 

Roughly  speaking,  both  the  total  airborne  dust  and 
dust  deposition  fall  on  the  same  approximate  curve.  This 
indicates  that  the  basic  aerodynamic  shear  dust  lofting  model 
can  be  characterized  with  a "universal"  curve.  There  is  a 
discrepancy  in  the  100  KT  results  between  10  and  100  secs, 
due  primarily  to  lower  speed  negative  phase  winds.  However, 
it  is  evident  that  one  will  not  lose  much  accuracy  if  a com- 
promise curve  is  used  for  this  time  period. 

The  rate  of  dust  deposition  is  comparable  to  the  dust 
lofting  rate  at  times  greater  than  100  secs.  Maximum  air- 
borne mass  can,  therefore,  be  estimated  from  these  curves  to 
be  0.4  MT  per  MT  yield. 

The  blast  wave  contribution  to  the  high  cloud  cannot 
be  represented  by  a single  curve,  as  evidenced  by  the  results 
of  the  three  calculations.  This  indicates  the  significance 
of  energy  per  unit  volume  dumped  into  the  fireball  and  the 
subsequent  influence  of  the  increasing  fireball  rise  velocity 
at  the  higher  yields.  One  could  anticipate,  therefore,  that 
cloud  structure  and  cloud  dynamics  are  closely  linked. 
Nevertheless,  it  would  appear  that  as  much  as  0.1  MT  of  dust 
per  MT  yield  in  the  high  cloud  can  be  traced  back  to  the  air 
blast  wave  contribution.  This  represents  approximately 
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Figure  3.9, 


Total  airborne  dust,  deposition  dust,  and  high 
cloud  dust  in  MT  dust/MT  yield  as  a function  of 
time  (scaled  to  1 MT)  for  three  SCOUR  surface 
burst  calculations  with  three  different  hydro 
flow  fields.  Note  that  these  results,  while 
qualitatively  consistent  with  subsequent  SCOUR2 
calculations,  are  superseded  by  results  pub- 
lished in  References  3 and  4. 


one-third  the  total  dust  believed  to  be  present  in  the  high 
altitude  cloud  produced  by  surface  bursts. 

3.3.2  Near-surface  Dust  Cloud 

The  extent  of  the  near-surface  dust  cloud  is  strongly 
dependent  on  the  yield  of  a surface  burst.  Figure  3.10  shows 
plots  of  dust  cloud  altitude  as  a function  of  range  from 
ground  zero  for  100  KT,  1 MT,  and  20  MT . In  addition,  the 
lesser  yield  cloud  altitude  versus  range  curves  are  scaled 
by  the  cube  root  of  the  yield  ratio,  (20/W)^^^,  for  direct 
comparison  to  the  20  MT  results.  Note  that  the  times  for 
which  the  dust  clouds  are  traced  are  equivalent  to  approxi- 
mately 120  secs  at  1 MT.  Cloud  altitude  versus  range  at 
180  secs  is  plotted  in  Figure  3.11  for  the  100  KT,  1 MT,  and 
20  MT  cases,  to  provide  a same-time  comparison  of  the  near- 
surface cloud  geometry. 

These  results  give  an  indication  that,  at  least  for 
the  early  part  of  the  negative  phase,  near-surface  cloud  de- 
velopment can  be  represented  by  the  time  history  of  a single  ; 

curve  that  bounds  the  three  results,  where  both  height  and  ! 

time  are  scaled  by  the  cube  root  of  the  yield  ratio.  ! 

The  potential  damage  to  be  caused  by  the  nuclear  dust  j 

storm  is  strongly  dependent  on  the  augmented  dynamic  pressure  i 

attributable  to  the  blowing  dust.  Figure  3.12  shows  the  peak  ; 

dust  dynamic  pressure  calculated  at  overpressure  loca-  | 

tions  up  to  100  psi.  Three  nuclear  bursts  are  compared  to  the  J 

MIXED  COMPANY  HE  calculation  (discussed  in  Chapter  IV) . Note 
that  the  nuclear  burst  calculations  are  only  mildly  influenced 
by  yield,  and  within  a factor  of  + 25  percent,  sur- 

face bursts  is  a function  only  of  overpressure  location. 

In  all  of  the  nuclear  cases,  the  peak  dust  dynamic  pres- 
sure was  achieved  by  dust  that  was  trailing  right  behind  the 
air  shock  at  an  altitude  of  1 meter  to  3 meters  above  the 
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the  spatial  dimensions  of  the  lower  yield  clouds  scaled  by  (20/W) -‘-/■i . 
Note  that  these  results,  while  qualitatively  consistent  with  subsequent 
SCOUR2  calculations,  are  superseded  by  results  published  in  References 
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Figure  3.12.  Peak  dust  dynamic  pressure  calculated  with 
SCOUR2  versus  peak  overpressure  for  three 
nuclear  surface  bursts  and  the  MIXED  COMPANY 
HE  detonation  (Chapter  IV) . Note  that  these 
results,  while  qualitatively  consistent  with 
subsequent  SCOUR2  calculations,  are  superseded 
by  results  published  in  References  3 and  4. 


ground  surface.  These  dust  particles  actually  preceded  the 
scouring  of  dust  off  the  surface  at  a given  radial  location. 

When  the  scale  of  the  SCOUR  calculation  is  reduced  to  the 
one-half  KT  HE  explosion,  the  high  density  dust  is  trapped 
in  the  boundary  layer  and  does  not  reach  the  free  stream 
velocities  necessary  to  achieve  the  higher  values  of 

The  total  peak  dynamic  pressure  loading  in  the  non- 
interactive formulation  of  the  nuclear  surface  burst  dust 

storm  is  the  sum  of  q,  . and  q . . These  values  are  plotted 

•"dust  ^air  ^ 

versus  overpressure  in  Figure  3.13.  They  are  compared  to  the 
light  dust  and  heavy  dust  total  dynamic  pressure  taken  from 
the  current  edition  of  DNA  Report  EM-1,  Part  I,  "Capabilities 
of  Nuclear  Weapons,"  1 July  1972.  The  basis  for  the  latter 
curves  is  obscure.  Presumably  they  are  developed  from  a 
limited  amount  of  dynamic  pressure  data  obtained  in  actual 
above-ground  nuclear  explosions.  In  most  cases,  these  were 
low  HOB  detonations,  and  DNA  Report  EM-1  dust  predictions  in- 
volve extrapolation  of  these  data  to  surface  bursts. 

SC0UR2  gives  lower  dynamic  pressure  predictions  at  the 
higher  overpressures,  possibly  due  to  the  absence  of  crater 
ejecta  in  current  SC0UR2  calculations.  Dynamic  pressures  for 
a subsequent  36.6  KT,  350  ft  HOB  calculation  were  also  com- 
pared to  DNA  Report  EM-1  values  (see  Reference  3),  and  these 
indicated  that  the  addition  of  neutron  lofted  dust  to  the 
region  just  under  the  blast  significantly  boosted  the  SC0UR2 
predictions  of  <3,pQrp*  The  higher  SCOUR2  dust  loading  and 
dynamic  pressures  at  the  lower  overpressures  (p  < 30  psi)  in 
Figure  3.13  may  be  attributable  to  lower  (ideal  surface)  air 
dynamic  pressures  used  in  the  derivation  of  the  DNA  curve. 

The  particle  size  distribution  in  the  near-surface  dust 
cloud  varies  with  time  and  altitude.  At  early  times  the  posi- 
tive phase  dust  cloud  closely  adheres  to  the  original  in  situ 
size  distribution.  However,  at  later  times,  the  large  particles 
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Total  Dynamic  Pressure,  ergs/cc 


Figure  3.13.  Peak  total  dynamic  pressure  versus  peak  over- 
pressure as  predicted  by  SC0UR2  and  DNA  EM-1 
for  nuclear  surface  burst.  Note  that  these 
results,  while  qualitatively  consistent  with 
subsequent  SC0UR2  calculations,  are  superseded 
by  results  published  in  References  3 and  4. 


i ' 

have  fallen  out  and  only  particles  with  diameters  less  than 
0.001  cm  are  left  in  the  near-surface  cloud.  Results  for  a 
1 MT  surface  burst  are  shown  in  Figure  3.14. 

Peak  dust  concentrations  in  the  positive  phase  range 
from  0. 1 X 10  ^ to  3 X 10  ^ g/cc  between  the  0.5  psi  and 
100  psi  overpressure  location.  In  the  negative  phase,  the 
more  diffuse  near-surface  cloud  consists  of  only  the  small 
size  particulates  (d^  < 0.001  cm)  in  concentrations  between  i 

1 X 10  * g/cc  and  2 x lo  g/cc.  ' 

The  total  mass  flux  of  positive  phase  dust  at  the  50  psi  i 

overpressure  station  is  plotted  versus  altitude  in  Figure  3.15. 

At  each  yield  there  is  an  altitude  (Y)  at  which  dust  activity 
is  most  intense.  Both  the  total  mass  flux  and  Y increased 
with  yield.  Note  that  for  the  one-half  KT  HE  calculation, 
intense  dust  activity  is  confined  to  altitudes  less  than 

2 meters  (Y  = 0.5  m) . At  multi-megaton  yields,  this  is  ex- 
panded  to  25  meters  (Y  = 7.5  m). 

-t+ 

The  peak  dust  impulse  in  the  positive  phase,  J p,u^dt, 

1/3  o a 

scaled  by  (W)  , is  plotted  as  a function  of  peak  overpressure 

in  Figure  3.16  to  illustrate  that  for  nuclear  surface  bursts, 

dust  impulse  (like  air  impulse)  can  be  scaled  by  the  cube  root 

of  the  yield  ratio.  This  is  because  the  peak  dust  dynamic 

pressure  is  almost  independent  of  yield  (see  Figure  3.12)  and 

1/3 

the  positive  phase  duration  scales  by  (W) 
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Airborne  Dust  Mass,  MT 


Figure 


3.14.  Airborne  dust  mass  for  each  particle  size  in  a 
SCOUR  calculation  of  the  air  blast  lofted  dust 
cloud  (1  MT  yield).  Note  that  these  results, 
while  qualitatively  consistent  with  subsequent 
SC0UR2  calculations,  are  superseded  by  results 
published  in  References  3 and  4. 
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Figure  3.15. 
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Altitude,  meters 

Positive  phase  dust  mass  flux  profiles  for  three 
nuclear  surface  burst  SCOUR  calculations  at  50  psi 
peak  overpressure  location.  Inset  gives  profile 
from  500  ton  HE  SCOUR  calculations  discussed  in 
Chapter  IV.  Note  that  these  results,  while  quali- 
tatively consistent  with  subsequent  SCOUR2  cal- 
culations, are  superseded  by  results  published  in 
References  3 and  4. 
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Figure  3.16.  Peak  dust  impulse  scaled  by  (W  in  MT)  as  a 

function  of  peak  overpressure  from  three  SCOUR 
calculations.  Note  that  these  results,  while 
qualitatively  consistent  with  subsequent  SCOUR2 
calculations,  are  superseded  by  results  published 
in  References  3 and  4. 


IV.  SCOUR  CALCULATIONS  OF  BLAST  WAVE 
DUST  LOFTING  EXPERIMENTS 

4.1  INTRODUCTION 

During  the  course  of  the  S^  effort  to  develop  the  SCOUR 
code,  two  major  atter..pts  were  made  to  obtain  blast  wave  dust 
storm  data  from  large  scale  high  explosive  (HE)  experimental 
detonations.  Dust  cloud  measurements  of  cloud  height  and 
particulate  concentration  as  a function  of  time  were  made  at 
MIDDLE  GUST,  EVENT  II (100  tons  of  TNT)  and  MIXED  COMPANY^^^^ 
(500  tons  of  TNT) . Although  the  data  recovered  were  incomplete 
and  far  from  comprehensive,  they  do  provide  an  experimental 
basis  for  comparJ.son  to  SCOUR  calculations. 

In  Reference  2,  the  previous  and  present  diffusion 
models  had  been  evaluated  by  comparing  the  results  of  calcula- 
tions achieved  with  a special  version  of  the  SCOUR  code  de- 
signed for  use  with  an  idealized  blast  wave  flow  field.  The 
focus  was  on  the  MIDDLE  GUST  results,  and  it  was  determined 
that  the  present  diffusion  model  was  adequate  for  the  computa- 
tion of  the  in  situ  dust  cloud  altitude  versus  time  as  well  as 
the  concentration  of  lofted  particulates.  (A  "soil"  of  5 mm 
diameter  glass  beads  was  used  in  the  MIDDLE  GUST  experiments 
to  determine  particulate  concentrations.)  A set  of  calcula- 
tions for  MIXED  COMPANY,  using  the  newer  diffusion  model 
(discussed  in  Reference  28)  were  made  shortly  thereafter,  and 
corroborated  some  preliminary  (unpublished)  results;  i.e., 
without  changing  the  SCOUR  parameters  from  those  used  in  the 
MIDDLE  GUST  calculations,  the  SCOUR  model  yielded  a high  degree 
of  correlation  to  the  experimental  measurements  of  the  in  situ 
dust  cloud. 

These  early  SCOUR  results  were  extremely  encouraging, 
but  due  to  the  incomplete  nature  of  the  experimental  results, 
it  was  only  possible  to  conclude  that  SCOUR  was  "in  the  right 
ballpark".  Moreover,  it  was  realized  at  the  end  of  the  present 
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report  period  that  a number  of  coding  errors  were  present 
in  the  special  version  of  the  SCOUR  code  used  in  these  cal- 
culations. While  the  errors  were  expected  to  have  a minor 
effect  on  overall  cloud  characteristics,  it  was  felt  that 
they  might  exert  a major  influence  on  the  SCOUR  predictions 
for  near-surface  particulate  concentrations,  particularly  in 
the  negative  phase.  The  corrections  to  the  SCOUR  code  were 
as  follows: 

• Boundary  layer  transition  — in  the  earlier  versions 
of  SCOUR,  transition  from  the  positive  phase  boundary 
layer  to  negative  phase  model  was  delayed  by  two  time 
cycles.  This  was  modified  so  that  the  negative  phase 
boundary  layer  calculation  is  initiated  immediately 
after  a negative  radial  (inward)  velocity  is  indi- 
cated for  a given  surface  zone. 

• Negative  phase  particle  lofting  — in  the  earliest 

version  of  SCOUR, particles  were  ejected  into  the 

flow  field  at  90®  to  the  horizontal.  A later  ver- 
r 2 1 

sion  included  a lofting  model  that  ejected  parti- 
cles, randomly,  between  10°  and  90°  to  th«=  horizontal. 
Due  to  a coding  error  (that  has  been  corrected) , the 
negative  phase  particles  were  ejected  in  opposition 
to  the  wind. 


• Random  number  sequence  — a coding  error  was  inadver- 

tently introduced  to  the  SCOUR  code  (1974  version) 
whereby  the  random  number  sequence  was  automatically 
restarted  w’henever  certain  computer  operations 
i occurred  during  a calculation.  This  results  in  a 

bias  in  the  random  numbers,  and  was  subsequently 

41 

corrected  for  the  SCOUR2  calculations.  ‘ ' 


Our  main  concern,  of  course,  was  to  determine  if  the 
tentative  conclusions  in  References  2 and  28  were  still  justi- 
fied; i.e.,  that  the  nominal  SCOUR  prescriptions  for  the  mass 
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lofting  coefficient,  = 0.1,  and  the  diffusion  parameter, 

k'  = 0.08,  were  adequate  for  the  calculation  of  near-surface 
particulate  concentrations  as  observed  in  the  MIDDLE  GUST  II 
(5  mn\  glass  beads)  and  MIXED  COMPANY  (in  situ  dust)  particu- 
late lofting  experiments. 

The  MIXED  COMPANY  event  was  recalculated  with  SCOUR2 
three  times  to  evaluate  the  effect  of  using  different  values 
for  the  negative  phase  diffusion  parameter.  The  results  were 
slightly  different  than  the  earlier  values  and  indicated  that 
the  nominal  prescription  was  superior  in  the  positive  and 
early  negative  phases  of  the  flow  field  (on  the  basis  of  com- 
parisons to  experiment) . 

Based  on  these  results  and  some  unanswered  questions 
about  the  earlier  MIDDLE  GUST  SCOUR  calculations,  the  cor- 
rected SCOUR2  model  was  then  used  to  recalculate  the  5 mm 
diameter  glass  beads  lofting  experiment.  These  new  results 
indicated  that  the  nominal  SCOUR  mass  lofting  coefficient 
(b^l^  = 0.1)  was  approximately  half  that  required  to  match 
the  experimental  results.  In  addition,  it  was  determined 
that  the  SCOUR  deposition  model  may  be  primarily  responsible 
for  the  overestimation  of  near-surface  particulate  concentra- 
tions at  late  times. 

In  this  chapter,  the  results  of  these  four  sets  of 
SCOUR  calculations  are  presented  in  the  following  sequence; 
MIDDLE  GUST  II  (1974),  MIXED  COMPANY  (1974),  MIXED  COMPANY 
(1975),  MIDDLE  GUST  II  (1975).  The  MIDDLE  GUST  II  (1974) 
results  have  been  published  in  Reference  2,  but  are  included 
in  this  chapter  for  ease  of  reference. 

4.2  MIDDLE  GUST  (1974  CALCULATIONS) 

The  MIDDLE  GUST  event  was  a detonation  of  a 100  ton 
sphere  of  TNT  located  two  sphere  radii  above  the  surface. 

An  assortment  of  particulate  lofting  experiments  were 


conducted  by  a TRW  research  team. These  consisted  of 
motion  picture  coverage  of  the  early  time  rate  of  rise  of 
the  in  situ  dust  cloud  kicked  up  by  the  air  shock  at  the 
50  psi  and  25  psi  overpressure  locations.  In  addition, 

5 mm  diameter  beads  were  placed  in  prepared  beds,  7 ft  by 
70  ft,  that  preceded  the  nominal  50  psi  and  25  psi  over- 
pressure stations.  Concentration  histories  of  5 mm  glass 
beads  at  the  25  psi  overpressure  location  were  measured  with 
a microwave  attenuation  system.  Concentration  data  were  not 
obtained  at  the  50  psi  location. 

The  SCOUR  dust  cloud  predictions  require  a hydro  flow 
field  input.  An  idealized  prescription  of  the  blast  wave 
flow  field  was  utilized,  to  avoid  complications  in  the  output 
of  hydro  calculations  of  HE  detonations  associated  v/ith  the 
presence  of  reaction  products.  (Moreover,  AFWL  tapes  of  the 
HE  calculations  made  available  to  do  not  cover  the  negative 
phase  of  the  event.) 


The  input  flow  field,  fully  described  in  Appendix  B, 

(331 

is  similar  to  that  used  by  Bannister.  Essentially,  the 

predicted  shock  location,  X^,  as  a function  of  time  t,  has 
been  fit  by  the  following  expression: 


(t) 


0.525 


(4.1) 


This  gives  the  shock  arrival  time  and  peak  (shock)  flow  con- 
ditions at  a given  radial  location,  X.  This  is  combined  with 
an  idealized  prescription  of  the  positive  phase  and  negative 
phase  flow  histories  at  X that  require  only  the  predicted 
positive  phase  duration  as  input.  The  air  is  assumed  to  be 
an  ideal  gas  (y  = 1.4)  and  detonation  products  are  not  in- 
cluded in  the  flow  field. 

Prior  to  the  detonation,  the  in  situ  soil  was  charac- 
terized as  a sandy  clay.  It  consisted  primarily  of  fine  grain 
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material.  The  measured  particle  size  distribution  is  given 
in  Figure  4.1.  For  SCOUR  calculations,  the  follov/ing  distri- 
bution was  used: 

Particle  Size  SCOUR  Characteristic  Mass  Fraction  of 

Range,  mm Particle  Diameter,  mm  Particle  Size  Class 


d 

< 0.030 

0.005 

40% 

0.030 

< d 

< 0.150 

0.005 

20% 

0.150 

< d 

< 0.400 

0.300 

20% 

0.400 

< d 

0.800 

20% 

4.2.1 

In 

Situ  Dust  Cloud 

The  peak  altitudes  of  particles  in  zones  above  the  50 
psi  and  25  psi  overpressure  locations  are  plotted  in  Figures 
4.2  and  4.3.  It  is  immediately  apparent  from  these  results 
that  there  is  a consistent  correlation  at  both  locations  be- 
tween the  observed  cloud  traces  from  the  two  cameras  and  the 
SCOUR  cloud  predictions.  The  smaller  diameter  particles 
(5  u,  70  u)  rise  at  a rate  intermediate  between  the  two  cloud 
rise  curves  taken  from  the  two  photographic  records.  At  both 
locations,  the  300  u particles  initially  rise  to  heights 
greater  than  the  observed  cloud.  Their  rate  of  rise  lessens, 
and  by  40-80  milliseconds  are  at  altitudes  comparable  to  the 
smaller  particles.  The  0.8  mm  particles,  delayed  by  the  time 
it  takes  to  achieve  lofting  velocities,  eventually  outclimb 
all  other  particle  sizes. 

The  concentration  profiles  at  specific  times  in  the 
calculation  are  plotted  (as  bar  graphs)  in  Figures  4.4a-c. 
These  are  compared  to  cloud  height  indications  from  the 
photographic  records.  It  is  seen  that  if  the  SCOUR  model  is 
accurate,  the  HYCAM  camera  was  blacked  out  by  concentrations 

- 4 - 4 

between  1 x 10  g/cc  and  3 x lo  g/cc  of  the  larger  grain 
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End  Positive  Phase 

Time,  milliseconds  (after  shock  arrival) 


Figure  4.2.  Peak  altitudes  of  SCOUR  (1974)  particles  above 
the  50  psi  overpressure  location  (IlIDDLE  GUST) 
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Concentration,  xlO'*  g/cc 

Figure  4.4.  Dust  concentration  levels  predicted  by  SCOUR  (1974) 
for  in  situ  MIDDLE  GUST  dust  cloud. 
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particles.  The  Millikan  camera  cloud  was  transparent  until 

— 4 

concentrations  of  small  diameter  particles  above  5 x 10  g/cc 
were  encountered. 

4.2.2  5 mm  Glass  Bead  Experiments 

A second  SCOUR  calculation  was  performed  where  the 
in  situ  soil  was  replaced  by  a new  soil  composed  of  5 mm  glass 
beads.  The  shear  transmission  coefficient  remained  set  at  0.1, 
and  no  other  model  modifications  were  made.  The  results  of 
this  calculation,  presented  in  Figures  4.5  and  4.6,  are  com- 
pared to  the  experimentally  obtained  concentration  histories 
at  the  6 inch  and  24  inch  altitudes  (25  psi  overpressure  loca- 
tion) . At  both  altitudes  there  is  excellent  agreement  between 
theory  and  experiment  for  peak  concentrations.  The  predicted 
concentration  time  history  at  the  6 inch  altitude  is  also  in 
good  agreement  with  experimental  results.  There  is  a signifi- 
cant difference  between  the  SCOUR  cloud  and  experiment  at  the 
24  inch  height  for  0.2  sec  < t < 1.0  sec,  perhaps  due  to  under- 
estimating the  diffusion  in  the  boundary  layer.  However,  the 
discrepancy  may  be  closely  related  to  uncertainties  in  simulat- 
ing the  late  time  flow  fields  as  well  as  the  aftereffects  of 
flow  asymmetries.  It  should  also  be  noted  that  at  late  times 
(t  > 1.0  secs),  the  SCOUR  concentrations  did  not  monotonically 


decay  as  the  experimental  saltation  values.  Rather,  the  SCOUR 
results  indicate  significant  particulate  saltation  and/or 
fallout  activity.  This  disagreement  is  most  likely  related  to 
the  particle  impact  deposition  model  (as  discussed  in  Section 
4.5)  . 


4.3 


MIXED  COMPANY  (1974) 


The  MIXED  COMPANY  HE  event  had  a spherical  charge  con- 
figuration of  500  tons  of  TNT,  centered  one  radius  above  ground 
zero.  Two  primary  measurements  were  made  to  study  the  dust 
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lofting  characteristics  of  the  air  blast  wave.  A concentra- 
tion time  history  of  the  in  situ  dust  cloud  was  measured  with 
laser  light  extinction  instrumentation  at  the  6 inch  height 
of  the  25  psi  overpressure  location.  Motion  picture  coverage 
of  the  early  time  dust  cloud  was  obtained  at  the  50  psi  over- 
pressure location.  In  addition,  supplemental  information  was 
acquired  at  the  25  psi  overpressure  position  from  post-test 
examination  of  7 foot  high  "greasy  stakes"  which  collected 
dust  impacting  on  the  four  flat  faces  of  the  stake.  These 
experiments,  performed  by  a TRW  research  team,  are  reported 
in  Reference  32. 

The  "dust"  at  the  MIXED  COMPANY  site  was  considerably 
modified  from  the  original  expectation  of  a silty  soil  similar 
to  the  MIDDLE  GUST  in  situ  dust.  Due  to  previous  snows,  the 
soil  was  quite  wet  and  cloddy  at  the  time  of  the  test.  Inter- 
pretation of  the  laser  light  extinction  data  required  that  the 
in  situ  particle  size  distribution  be  deduced  from  the  dust 
captured  by  the  greasy  stake.  This  distribution,  indicated 
by  the  x's  in  Figure  4.7,  v^as  simplified  into  a four  particle 
size  in  situ  soil  for  the  SCOUR  calculation  of  the  MIXED 
COMPANY  dust  cloud. 


Particle  Diameter 
0 . 03  mm 
0.20  mm 

1 . 0 mm 

3 . 0 mm 


Mass  Fraction 
30% 

30% 

20% 

20% 


An  idealized  flow  field  was  used  to  simulate  the  air 
blast  wave  passage  and  subsequent  flow.  It  is  a scaled  ver- 
sion of  the  MIDDLE  GUST  flow  field.  Surface  zones  were  4.275 
meters  wide  for  this  calculation  (scaled  from  MIDDLE  GUST  by 
the  cube  root  of  the  yield  ratio,  see  Appendix  B) . 


i 
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Figure  4.7.  Particle  sizes  used  in  SCOUR  simulation  of 

MIXED  COMPANY  event  compared  to  distribution 
used  to  reduce  laser  light  extinction  data [32]. 
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4.3.1  In  Situ  Dust  Concentration 


SCOUR  predictions  of  the  dust  concentration  at  the 
6 inch  height  are  compared  with  experimental  results  in 
Figure  4.8.  Values  at  the  35  psi  and  25  psi  overpressure 
location  are  given  to  provide  an  indication  of  the  spread 
in  concentration  data  predicted  by  the  SCOUR  code  for  this 
HE  event. 

In  general,  the  predictions  with  the  standard  SCOUR 
parameters  are  in  good  agreement  with  the  experimental  re- 
sults (+  1 X 10**  g/cc)  . However,  having  data  at  only  one 
altitude  and  at  one  overpressure  location,  does  not  provide 
adequate  validation  of  the  SCOUR  model.  The  next  figure 
(4.9)  shows  the  SCOUR  predictions  for  different  values  of 
the  diffusion  velocity  parameter  (k^  = 0.02,  0.08,  and  0.16). 
It  is  evident  that  the  nominal  prescription  (k^  = 0.08) 
matches  the  data  best  in  the  late  time  regime  (t  > 0.5  sec) 
but  may  be  over-diffusing  dust  at  early  times. 

To  more  appropriately  validate  the  SCOUR  model,  dust 
profile  information  is  required.  Figure  4.10  compares  the 
predicted  dust  concentrations  between  35  cm  and  70  cm  above 
the  original  surface  as  a function  of  time  for  cases  computed 
with  different  values  of  k' . Agreement  with  experiment  in 
this  altitudinal  range  (if  data  were  available)  would  greatly 
strengthen  the  argument  for  a particular  set  of  k'  values. 

It  would  also  provide  greater  insight  into  the  applicability 
of  the  aerodynamic  mass  lofting  prescription. 

4.3.2  Dust  Cloud  Height 

The  dust  cloud  height  versus  time,  measured  with  photo- 
graphic techniques  at  the  50  psi  overpressure  location,  is 
compared  to  the  maximum  particle  altitudes  predicted  by  SCOUR 
in  Figure  4.11.  The  various  symbols  represent  the  different 
particle  sizes  used  in  the  SCOUR  calculation. 
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Figure  4.8.  Comparison  of  calculated  SCOUR  (1974)  dust  con- 
centrations at  the  35  psi  and  25  psi  peak  over- 
pressure locations  to  MIXED  COMPANY  measure- 
ments [32]  (at  nominal  25  psi  position). 
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Figure  4.9.  Comparison  of  SCOUR  (1974)  dust  concentrations 
at  the  25  psi  overpressure  location  for  three 
values  of  diffusion  coefficient,  k',  and  the 
experimental  observations  made  during  MIXED 
COMPANY  event. [32] 
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Figure  4.10.  SCOUR  (1974)  predictions  of  dust  concentration 

time  histories  between  35  cm  and  70  cm  above  the 
ground  surface  at  the  25  psi  overpressure  loca- 
tion for  a 500  ton  HE  explosion,  using  three 
different  diffusion  coefficients. 
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Figure  4.11.  Maximum  height  of  various  size  particles  in 

SCOUR  (1974)  calculated  dust  cloud  compared  to 
photographic  observations,  MIXED  COMPANY. [32] 
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From  this  comparison  it  is  evident  that  the  smaller 
particles  (d^  = 0.03  mm,  0.20  mm)  are  embedded  in  the  early 
time  dust  cloud  while  the  larger  particles  (1  mm,  3 mm)  indi- 
cate the  outer  periphery  of  the  dust.  Significantly,  the 
standard  k'  prescription  diffuses  particles  out  to  altitudes 
of  one  meter  by  the  end  of  the  positive  phase,  while  reduced 
positive  phase  diffusion  (k'’  = 0.02)  limits  the  small  particle 

altitudes  to  about  0.5  meter  (plus  a slower  initial  cloud 

r 321 

rise).  The  greasy  stakes  experiment  indicated  that  the 

small  diameter  particles  were  concentrated  throughout  the 
bottom  meter,  giving  supportive  evidence  for  the  higher  dif- 
fusional  rates  behind  the  blast  wave. 

4.4  MIXED  COMPANY  CALCULATIONS  (1975) 

The  1975  SCOUR  calculations  of  the  MIXED  COMPANY  event 
were  made  with  the  same  idealized  flow  field  used  in  the 
earlier  calculations.  In  each  of  the  three  new  calculations, 
the  mass  lofting  parameter  was  fixed  at  0.1  and  the  positive 
phase  diffusion  parameter  was  kp  = 0.08.  The  only  change  in 
each  case  was  the  negative  phase  diffusion  parameter;  i.e. , 
k^  = 0.02,  0.08,  0.16.  (The  kp  = 0.08  combination  is  the 
nominal  prescription.) 

Dust  concentrations  at  the  25  psi  overpressure  location 
in  each  of  these  calculations  are  plotted  versus  time  in  Fig- 
ure 4.12.  They  have  been  computed  by  averaging  the  dust  con- 
centration between  10  cm  and  35  cm  above  the  surface  for  three 
radial  zones  (213.7  m ^ r £ 226.6  m)  . Also  plotted  in  Figure 
4.12  are  the  experimental  measurements  of  dust  concentration 
at  a 6 inch  altitude  reported  in  Reference  32. 

It  is  clear  from  the  comparison  of  SCOUR  predictions  to 
the  experimental  results  that  SCOUR  calculates  dust  concen- 
trations that  are  in  qualitative  agreement  with  experiment  for 
any  value  of  k^.  There  are  differences  between  the  1975 
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Dust  concentrations  between  10  cm  and  35  cm 
above  the  original  surface  at  the  25  psi  over 
pressure  location  of  the  MIXED  COMPANY  event, 
computed  with  kp  = 0.08  and  three  values  of 
(SCOUR  1975).  Results  are  compared  to  ex- 
perimental observations  made  at  15.2  cm  alti- 


nominal  prescription  results  and  those  presented  in  Figure 
4.9  (1974),  most  of  which  are  inconsequential.  The  major 
change  in  the  1975  results  is  in  the  predictions  for  the 
late  negative  phase.  It  is  seen  in  Figure  4.12  that  the 
nominal  prescription  yields  an  increase  in  concentration  as 
the  wind  speed  decays.  This  is  contrasted  by  the  1974  re- 
sults which  were  stabilized  in  a range  of  concentrations  that 
very  nearly  matched  the  experimental  values.  Moreover,  the 
other  values  of  negative  phase  diffusion  used  in  the  1975 
calculations,  i.e.,  k'  = 0.02,  0.16  provided  better  agree- 
ment with  experiment  during  the  late  negative  phase  0.6  sec 
< t < 1.0  secs. 

It  should  also  be  noted  that  for  the  positive-early 
negative  phase  time  period,  0.2  sec  < t < 0.6  sec,  the  nominal 
prescription  in  the  1975  calculations  provided  the  closest 
agreement  with  experimental  results,  and  represented  an  im- 
provement over  the  earlier  calculation  (see  Figure  4.9). 

On  the  basis  of  these  calculations,  the  nominal  pre- 
scription was  retained  for  the  final  set  of  SCOUR  surface 
burst  calculations,  the  results  of  which  are  presented  in 
References  3 and  4.  It  was  felt  that  better  positive  phase 
agreement  was  essential  because  it  is  the  positive-early 
negative  phase  dust  that  contributes  to  the  most  intense 
period  of  the  nuclear  dust  storm  and  is  the  primary  source 
of  dust  for  the  late  time  dust  cloud. 

However,  the  discrepancy  between  prediction  and  ex- 
periment at  late  times  in  the  negative  phase  is  a matter  of 
concern  because  it  indicates  the  possibility  that  the  SCOUR 
model  may  be  inadequate  during  this  period.  Here  it  should 
be  stated  that  the  ideal  flow  used  in  the  SCOUR  calculations 
does  not  include  the  effects  of  blast  wave  asymmetries, 
natural  afterwinds,  or  excess  turbulence  in  the  decay  process. 
All  of  these  could  play  an  important  role  in  the  suspension 
capacity  of  the  late  time  dust  flow  field. 
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Another  point  that  had  not  been  previously  considered 
was  that  the  SCOUR  model  includes  a dust  deposition  model 
that  has  not  been  validated.  In  SCOUR,  particles  that  impact 
the  surface  are  assumed  to  eject  other  particles  of  equal 
mass  with  roughly  10  percent  of  the  impact  momentum.  If  the 
ejected  particle  velocity  would  be  less  than  that  which  is 
required  for  a particle  to  maintain  saltation  (as  noted  in 
Chapter  II,  see  Figure  2.2),  a new  particle  is  not  ejected 
and  deposition  occurs  at  the  point  of  impact. 

Examination  of  the  impact  data  in  the  1975  calcula- 
tions reveals  that  deposition  begins  at  about  t = 0.75  secs 
after  detonation.  Thus,  in  the  late  negative  phase  we  could 
anticipate  that  the  SCOUR  deposition  model  plays  a major 
role  in  the  near-surface  dust  concentration  predictions. 

This  is  corroborated  in  the  average  concentration  (at  17.5 
cm)  versus  time  plots  for  the  three  calculations  in  Figure 
4.13.  These  differ  from  those  presented  in  Figure  4.12  in 
that  the  concentrations  from  the  bottom  10  cm  have  been 
averaged  with  those  in  the  altitude  range  10  cm  to  35  m. 

The  mean  height  (17.5  cm)  is  very  close  to  the  15.2  cm  alti- 
tude at  which  the  experimental  observations  were  made.  After 
t = 0.55  second,  the  concentrations  rise  steadily,  indicating 
the  strong  possibility  that  saltating  particles  in  the  near- 
surface region  could  be  strongly  influencing  the  near-surface 
dust  concentrations.  This  may  be  an  especially  important  ef- 
fect in  the  MIXED  COMPANY  dust  lofting  results  because  the 
in  situ  soil  was  damp  and  cloddy.  Thus,  it  is  possible  that 
late  negative  phase  saltation  was  impeded  by  these  surface 
conditions. 

4.5  MIDDLE  GUST  II  (1975) 

The  1974  MIDDLE  GUST  SCOUR  calculations,  presented  in 
Figures  4,2,  4.3,  4.4,  and  4.5,  indicated  that  the  revised 


Time  (since  detonation),  sec 

Figure  4.13.  Dust  concentrations  between  0 and  35  cm  above 

the  original  surface  at  the  25  psi  overpressure 
location  of  the  MIXED  COMPANY  event,  computed 
with  kp  = 0.08  and  three  different  values  of 
(SCOUR  1975) . Results  are  compared  to  experi- 
mental observations  at  15.2  cm  altitude. 


SCOUR  models  for  particulate  lofting  and  diffusion  (incor- 
porated in  SCOUR  in  1973)  provided  an  adequate  correlation 
with  the  near-surface  observations  of  the  blast  wave  dust 
cloud  rate  of  rise  and  the  concentrations  of  5 mm  diameter 
glass  beads  that  were  lofted  by  the  windstorm  at  the  25  psi 
overpressure  location.  However,  the  code  deficiencies  cited 
in  Section  4.1  necessitated  a recalculation  of  the  event. 

In  addition,  the  discrepancy  between  the  late  time  SCOUR  pre- 
dictions of  high  bead  concentrations  (for  t > 2 secs)  and  the 
near  zero  concentration  readings  of  the  experimental  apparatus 
remained  to  be  better  understood.  A second  area  of  concern 
was  that  the  problem  setup  for  the  1974  calculations  assumed 
a continuous  bead  surface  that  extended  beyond  the  actual  ob- 
servation point  used  for  comparison  to  experiment.  This  was 
in  slight  contrast  to  the  experimental  setup  where  bead  con- 
centrations were  actually  observed  slightly  dovmstream  of  the 
end  of  the  bed. 


Four  SC0UR2  calculations  were  made  to  determine  if  the 
above-mentioned  problem  areas  were  important  to  the  quantita- 
tive and  qualitative  nature  of  the  SC0UR2  bead  concentration 
predictions  for  the  MIDDLE  GUST  II  event.  The  particle  ejec- 
tion routine  for  the  negative  phase  was  corrected  and  particle 
ejection  was  limited  to  the  radial  extent  of  the  two  glass  bead 
beds  used  in  the  experiment.  Bead  concentration  was  monitored 
at  the  surface  zones  near  the  end  of  the  bed  in  three  vertical 
bins . 


Vertical  Bin 
1 
2 
3 


Altitude  Range 
0 — 10  cm 
10  - 35  cm 
35  — 70  cm 


Mean  Altitude 
5.0  cm 

22.5  cm 

52.5  cm 


Each  surface  zone  covered  2.5  meters  in  the  radial  direction. 
As  observed  in  previous  SCOUR  calculations,  it  is  necessary 
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to  average  over  three  radial  zones  to  obtain  relatively  smooth 
curves  for  concentration.  For  the  present  calculations,  SCOUR 
predictions  were  made  by  averaging  over  the  last  two  zones  in 
the  glass  bead  bed  and  the  first  downstream  "inactive"  zone 
(radial  zones  44-46)  . 

The  results  of  the  first  1975  calculation  were  not 
nearly  as  satisfactory  as  those  indicated  in  Figures  4.5  and 
4.6.  Concentrations  in  the  bottom  three  vertical  bins  are 
compared  in  Figure  4.14  to  the  experimental  measurements  made 
at  heights  of  15  cm  and  60  cm.  The  concentrations  are  approxi- 
mately one-third  the  observed  values,  and  do  not  decay  signifi- 
cantly at  late  times.  A second  set  of  concentrations  from  the 
same  calculation  are  plotted  in  Figure  4.15.  In  this  case, 
the  concentrations  were  computed  by  averaging  over  the  last 
three  radial  zones  of  the  bed  (all  "active"  bead  lofting  zones) . 
The  concentrations  are  significantly  higher,  but  still  below 
the  experimental  values. 

To  improve  SC0UR2  predictions  for  the  "beads"  experi- 
ments, the  angle  for  particle  ejection  in  the  SCOUR  mass 
lofting  routine  was  limited  to  angles  between  10°  and  40°. 

The  results  for  this  second  calculation  are  shown  in  Figure 
4.16.  This  had  the  beneficial  effect  of  raising  bead  concen- 
trations without  qualitatively  changing  the  temporal  history. 
Note  that  at  the  lower  altitudes,  the  apparent  (first)  peak  in 
concentration  is  predicted  to  occur  at  about  0.9  sec  (0.78  sec 
after  shock  arrival).  For  the  third  vertical  zone,  the  first 
peak  occurs  at  about  0.5  sec  (0.38  sec  after  shock  arrival). 

Both  of  these  are  in  excellent  agreement  with  peak  concentra- 
tion times  observed  during  MIDDLE  GUST.  A discrepancy  in 
concentration  magnitude  still  remains,  however,  and  the  mono- 
tonic decay  of  the  experimental  data  beyond  these  times  is 
not  simulated. 
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Results  of  first  SC0UR2  (1975)  calculation 
(using  nominal  prescription)  compared  to 
MIDDLE  GUST  II  experimental  results.  [31] 


77 


middle  gust  II  SQOUR  Concentxations 
25  psi  Overpressure  Averaged  Over  Radial 
location  Zones  43-45. 

SCOUR  Bin  Altitudes 


MeasuTf 
Values, 
24  in. 
Alti-  . 
tude  ' 


Resu.lts  of  first  SCOUR2  (1975)  calculation 
using  nominal  parametric  prescription  but 
averaging  over  last  three  "active"  surface 
zones,  compared  to  MIDDLE  GUST  II  experi- 
mental results. [31] 


Figure  4.16.  Results  of  second  SC0UR2  (1975)  calculation. 

The  range  of  particle  ejection  angles  was  re 
duced  to  10°  < 0 < 40°. 


since  we  are  interested  in  concentrations  at  a 15  cm 
(6  inch)  height,  and  the  gradients  in  this  region  are  quite 
high  (as  indicated  by  the  large  difference  between  concentra- 
tions in  vertical  bins  1 and  2) , the  average  dust  concentra- 
tion throughout  zones  1 and  2 (0  £ y ^ 35  cm,  y = 17.5  cm)  is 
a more  appropriate  comparison  to  the  experimental  data.  This 
average  for  the  "tight  angle"  calculation  is  shown  in  Figure 
4.17. 

The  periodicity  in  the  concentration  time  histories 
(most  evident  in  Figure  4.17)  was  thought  to  be  associated 
with  the  bouncing  of  particles.  Inspection  of  particle  impact 
data  corroborated  this  suspicion.  The  515  SCOUR2  particles 
that  were  lofted  by  the  air  shock  underwent  over  1500  impacts 
in  three  seconds.  To  determine  the  effect  of  eliminating 
particle  bouncing  from  the  SCOUR2  predictions,  the  minimum  im- 
pact velocity  necessary  to  eject  a surface  particle  by  a de- 
scending particle  was  raised  by  a factor  of  five  over  the 
nominal  SCOUR  prescription.  This  had  the  result  of  "killing" 
three  out  of  every  four  SCOUR2  particles  that  impacted  the 
ground.  ("Killed"  particles  are  deposited  at  the  point  of 
impact.)  The  concentration  time  histories  for  the  three 
vertical  zones  with  the  impact  modifications  in  this  third 
SCOUR2  calculation  are  compared  to  the  experimental  results 
in  Figures  4.18  and  4.19.  The  decrease  in  particle  bouncing 
led  to  predictions  of  time  to  peak  concentration  that  are  in 
agreement  with  the  experimental  observations.  However,  the 
magnitude  of  glass  bead  particle  concentrations  predicted  by 
SCOUR2  are  lower  by  a factor  of  three,  and  the  decay  to  zero 
concentration  is  too  rapid. 

In  the  subsequent  fourth  SCOUR2  calculation,  the  momen- 
tum transmission  coefficient  between  impacting  and  (to  be) 
ejected  particles  was  increased  from  0.1  to  0.7.  This  results 
in  more  particle  impact  lofting  as  well  as  a more  energetic 
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Results  of  the  third  SCOUR2  (1975)  calculation, 
y = 52.5  cm,  using  10°  ^ 6 ;!  40°  and  multiplying 
the  minimum  impact  velocity  to  loft  an  additional 
particle  by  a factor  of  five. 


saltation  layer  near  the  surface.  The  SC0UR2  concentrations 
at  y = 17.5  cm  and  y = 52.5  cm  are  compared  to  experimental 
results  in  Figures  4.20  and  4.21.  The  SCOUR2  concentrations 
averaged  over  the  three  "live"  bed  zones  (zones  43-45)  are 
also  plotted  for  comparison  to  the  nominal  case  (zones  44-46) 

I (the  last  two  "live"  zones,  and  the  first  "dead"  zone  just 

downstream  of  the  bed) . It  is  clear  that  the  SCOUR2  concen- 
trations are  still  lower  than  experiment,  but  the  time 
I variations  are  similar  to  the  experimental  data.  If  the  mass 

f lofting  coefficient  is  multiplied  by  a factor  of  2.3  (to 

(roughly  match  peak  concentrations  with  experiment)  the  SCOUR2 

predictions  are  in  excellent  agreement  with  the  MIDDLE  GUST 
I data.  This  is  demonstrated  in  Figures  4.22  and  4.23. 

If  we  accept  the  experimental  data  for  MIDDLE  GUST  as 
accurate,  the  results  of  these  parametric  studies  with  the 
SCOUR2  code  are  indicative  that  the  glass  beads  surface  medium 
was  more  susceptible  to  aerodynamic  lofting  than  actual  soil. 
In  addition,  as  indicated  by  shock  tube  studies,  our  calcula- 
tions show  that  the  glass  bead  particles  may  be  limited  to 
ejection  angles  of  less  than  40°.  It  is  clear  that  particle 
impact  conditions  strongly  influence  the  bead  "cloud"  dynamics 
and  need  to  be  more  carefully  studied.  An  encouraging  result 
is  that  the  basic  time  variation  of  the  glass  beads  concen- 
‘ tration  is  predicted  by  the  SCOUR2  calculations.  Hence,  it 

appears  that  the  dynamics  of  "bead"  lofting  processes  can  be 
well  represented  by  the  SCOUR2  code  models  and  numerical  pro- 
cedures. 


Figure  4.22.  Results  of  the  fourth  SC0UR2  (1975)  calculation,  y = 17.5  cm,  modified 
by  increasing  the  aerodynamic  lofting  coefficient  from  0.1  to  0.23. 
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V.  FLUIDIZED  BED  CALCULATIONS  OF  THE 

REVERSE  PERCOLATION  PROCESS 

5.1  INTRODUCTION  AND  SUMMARY 

It  was  recognized  early  in  dust  lofting  research  that 
a possible  mechanism  for  dust  blowoff  from  the  surface  was 
the  so-called  reverse  percolation  process. ^ ^ The  passage 
of  a high  pressure,  air  blast-wave  pulse  over  the  surface  of 
a porous  soil  could  lead  to  the  pressurization  of  air  trapped 
in  the  pores.  As  the  overpressure  decreases  above  the  soil, 
the  (high  pressure)  trapped  air  would  flow  upwards  through 
the  soil  matrix  (i.e.,  reverse  percolation),  entraining  soil 
particles  if  there  was  a strong  enough  pressure  gradient. 

Reverse  percolation  was  observed  in  a series  of  simple 

r 8 1 

experiments  performed  by  Shock  Hydrodynamics  in  1973.^  ‘ In 
these  tests,  a 20  cm  diameter  column  of  packed  sand  (650  ym 
diameter  (mean) , 33  percent  porous)  was  exposed  to  surface 
overpressures  that  corresponded  to  the  100  psi  overpressure 
environment  for  a 1 MT  surface  burst.  Significant  fluidiza- 
tion and  upward  motion  of  the  near-surface  sand  was  observed, 
and  a large  fraction  of  the  180  cm  long  column  was  set  in 
motion. 

r 2 1 

In  an  earlier  report,  ^ s’  performed  a simplified 
analysis  of  the  reverse  percolation  process  at  low  overpres- 
sures and  determined  that  it  was  possible  to  fluidize  the 
surface  particulates,  but  the  flow  velocity  was  too  small  to 
actually  eject  particles  into  the  positive  phase  air  flow. 

It  was  believed,  hov;ever,  that  additional  research  should  be 

performed  on  the  degree  to  which  near-surface,  porous  layers 

could  be  fluidized  and/or  ejected  into  the  negative  phase  air  ; 

flow. 

Relevant  to  the  dust  lofting  problem,  had  initiated 
the  research  and  development  of  a new  computer  code,  FLUB 
(FLUidized  Bed  code) . This  work  was  originally  performed  in 
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S 's  IR&D  program.  ' The  intent  was  to  better  under- 

stand the  dynamics  of  fluidized  beds,  and  to  develop  ad- 
vanced numerical  codes  that  could  be  used  to  study  fluidized 
bed  problems.  The  FLUB  code,  described  in  the  next  section 
(5.2)  has  been  used  to  calculate  the  fluidized  bed  dynamics 
observed  in  the  shock  tube  experiments ^ and  to  make  some 
preliminary  estimates  of  the  porous  soil  response  at  a high 
overpressure  position  under  nuclear  surface  burst  conditions. 

The  FLUB  results,  described  in  Sections  5.3  and  5.4, 
demonstrate  that  the  air  trapped  in  pores  of  the  soil  matrix 
will  fluidize  porous,  strengthless  soils  down  to  depths  of 
tens  of  centimeters.  Large  quantities  of  dust  can  be  given 
positive  vertical  velocities.  If  the  simplified  model  and 
assumptions  that  underlie  the  application  of  the  FLUB  code 
are  valid,  large  quantities  of  dust  will  be  given  small  (but 
finite)  positive  vertical  velocities.  Thus,  the  reverse  per- 
colation effect  could  dramatically  alter  the  air  flow/surface 
interface  conditions  from  that  generally  assumed  for  aero- 
dynamic lofting.  Namely,  the  low  lying  saltation  layer  would 
be  replaced  by  a fluidized  bed  with  much  higher  particulate 
concentrations.  Presumably,  the  mass  entrainment  at  the  air 
flow/fluidized  soil  interface  would  still  be  proportional  to 
the  shear  at  the  interface,  and  the  shear  may  still  be  modeled 
by  a rough  wall  shear  law  (as  in  SC0UR2) . Thus,  the  simple 
models  for  mass  lofting  and  shear  presently  in  SC0UR2  may  be 
an  appropriate  formulation/conceptualization  for  entrainment 
from  fluidized  beds,  but  the  parameters  may  require  reevalua- 
tion. 

Although  the  FLUB  calculations  provide  significant  in- 
sight to  a potential  problem  area  in  the  SCOUR  code  formula- 
tion, this  line  of  study  and  model  development  was  not  further 
explored  in  the  study.  The  chief  reason  was  that  a thorough 
understanding  of  the  reverse  percolation  process  would  involve 


90 


a far  more  definitive  treatment  of  soil  models  than  that 
possible  within  the  scope  of  the  SC0UR2  development  program. 
Such  an  effort  would  include  porosity,  saturation,  strength 
effects,  and  soil  matrix  analyses.  Ultimately,  this  would 
involve  parametric  code  calculations  to  gauge  model  sensiti- 
vity and  additional  experimentation  to  validate  model  "num- 
bers" . 

Secondly,  the  reverse  percolation  process  is  signifi- 
cant only  in  the  negative  phase  of  the  blast  wave  flow  field, 
over  a limited  range  and  time.  The  SC0UR2  model  for  aero- 
dynamic scouring  by  the  positive  phase  blast  wave  would  not 
be  significantly  affected,  hence  positive  phase  dust  loads 
on  structures  would  not  be  affected.  Taking  this  into  ac- 
count, and  realizing  that  a soil  matrix  with  small  but  finite 
cohesion  would  seriously  reduce  the  reverse  percolation  loft- 
ing potential,  this  topic  was  given  a lower  priority  than  the 
popcorning  due  to  thermal  radiation  and  neutron  absorption  by 
groundwater. 

5 . 2 FLUB 

The  FLUB  code  was  originally  developed  to  solve  the 
governing  equations  for  an  unsteady,  chemically  inert  gas- 
fluidized  bed  in  one  space  dimension.  These  equations,  de- 
veloped within  the  framework  of  the  Theory  of  Interacting 
Continua  (TING) , express  the  principles  of  conservation  of 
mass  and  momentum  for  both  components  (solid  particles  and 
interstitial  gas)  of  the  mixture.  Derivations  of  these  equa- 
tions can  be  found  in  Garg. is  in  the  process  of  ex- 
tending the  methodology  to  two  space  dimensions,  including 
chemical  reactions,  energy  transfer,  and  multiple  species. 

In  the  original  version  of  FLUB,  the  fluid  is  treated 
as  an  ideal  gas,  which  may  be  either  isothermal  or  isentropic. 
The  solid  is  considered  to  consist  of  a large  number  of 
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uniform-size  spherical  particles;  this  collection  of 
particles  behaves  like  a Newtonian  fluid  with  bulk  and 
shear  "viscosities"  which  may  depend  in  an  arbitrary  way 
upon  the  local  solid  concentration.  Each  individual 
particle  is  treated  as  an  incompressible  sphere,  but,  in 
the  aggregate,  the  particles  may  be  pressurized  if  the 
local  packing  density  exceeds  a preset  value.  The  solid 
and  fluid  phases  are  coupled  through  a friction  law  which 
approaches  Stokes  law  in  the  thin-mixture  limit  and  which 
approximates  flow  through  a porous  matrix  at  high  packing 
densities.  The  effective  "permeability"  as  a function  of 
porosity  was  matched  to  data  collected  by  Jackson in 
the  porosity  range  0.4-0. 5,  as  well  as  to  Stokes  law. 

The  numerical  procedure  employs  rather  elaborate 
implicit  iterative  techniques  to  solve  the  highly  nonlinear 
(and  occasionally  singular)  system  of  equations  which  re- 
sults from  using  realistic  constitutive  relations  in  the 
fundamental  governing  equations.  The  scheme  is  basically 
a finite  difference  procedure  in  both  the  space  and  time 
domains.  Simultaneous  use  of  an  Eulerian  and  a Lagrangian 
grid  permits  free  surfaces  to  be  accurately  delineated. 

The  Lagrangian  features  in  no  way  depend  on  the  ordering  of 
the  Lagrangian  nodes,  so  that  the  scheme  may  be  readily 
generalized  to  multidimensional  geometries  without  restric- 
tions upon  grid  distortion. 

5.2.1  Fundamental  Relations 

r 34  ] 

The  field  equations  developed  by  Garg  for  the  one- 

dimensional case  may  be  expressed  as  follows: 

1.  Conservation  of  Solid  Mass: 

II  (6u)  * 0 (5.1) 
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2. 


Conservation  of  Gas  Mass: 


It  + fe  [(i-e)PjUji  = 0 


Conservation  of  Solid  Momentum: 


(5.2) 


It  IjT  ■ 

- h '«> 

1 

" Ps 

3P 

3x 

i 

1 

[ ^ 

(5.3)  j 

9x 

Ps 

3x 

4. 


Here, 


Conservation  of  Gas  Momentum: 
k(e)  3P 

UjT  = U - -jz. — 

f (i-6)nf  3x 

pg  = solid  density  (a  constant) 
0 = 1-  ())  = 1-  porosity 
Pj  = gas  density 
= solid  velocity 


(5.4) 


u^  = gas  velocity 

A.,y  = "viscosity"  coefficients  for  solid  (depend  on  6) 

= gas  viscosity 

g = gravity  acceleration 

P = gas  pressure  (interstitial) 

0f  = solid  kinematic  pressure 

k = solid  permeability  (depends  on  6 and  the  nature 
of  the  particulate) 


The  gas  is  taken  as  polytropic,  so: 
P = (y-DC^T  p^ 


(5.5) 
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v’here 


= specific  heat  at  constant  volume 


T = gas  temperature 

Y = ratio  of  specific  heats  for  gas. 

Tv/o  forms  of  gas  thermodynamics  will  be  considered; 
Isothermal,  in  which  case: 


T = T = const 


(5.6a) 


so  that  the  equation  of  state  becomes 


p = (y-DV^p^ 


and 


Adiabatic,  in  which  case: 


(5.7a) 


(5.6b) 


SO  that 


\T-1 


P = (Y-l)C, 


7 0 V Pf Q / 


(5.7b) 


Here,  T and  are  the  reference  state  designating  the 

0 to 

particular  isentrope  under  consideration. 

Tlie  gas  viscosity  is  taken  as  proportional  to  absolute 
temperature,  so 


Hf  = nfo(T/T^) 


Thus,  for  the  isothermal  case. 


'^f  ’^fO 
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AO>A030  114 

unclassified 


systems  science  and  software  la  JOLLA  CALiF 
NEAR'SURFACE  nuclear  oust  cloud  studies. (U) 

JAN  77  J W KIRSCH  ONA001-73-C-0128 

SSS'R-76-2826  DnA>4332F  NL 


and  for  the  adiabatic  case 

hf  = nfo(Pf/Pfo)^"^  (5.8b) 

The  kinematic  solid  pressure  6f  is  expressed  in  the  following 
mathematical  form: 

0f  = ^ (6-6  ) 2 for  6 > 6 

2 0 0 (5.9) 

= 0 otherwise 

where  a and  0 are  specified  constants. 

0 

The  ratio  of  the  permeability  (k)  to  the  square  of  the 
solid  particle  radius  is  simply  a function  of  6: 

K(e)  = k/r^  . (5.10) 

0 

The  following  form  has  been  found  to  fit  available  data  for 

f 3 6 1 

spherical  particles  in  the  y^^^nge  0.4  to  0.5,  and  also 

matches  Stokes  law  in  the  limit  as  0 approaches  zero: 


K(0) 


2 

9 0 


X exp 


(-4.093  + O.3625[l-0]) 


(5.11) 


(see  Figure  5.1).  For  the  solid  viscosity  terms,  it  is  con- 
venient to  define: 


v(0) 


(5.12) 


and  to  use  the  form: 


(5.13) 

the  kinematic  viscosity  measured 
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k(6) 


R Stokes  Law  Asymptote 


Present  Fit 
2 1-0 


/ ic  = I exp  j [-4.093+0.365 
(1-0)] 


Measureroents 
(Jackson, 
Reference  8] 


0.001 


0.1  0.2  0.3 


0.4  0.5  0. 


Fiqure  5.1.  Dimensionless  permeability  {<  - k/r^)  as  a 

function  of  porosity  (=1-0)  fitted  to  Stokes 
Law  and  Jackson  [36]  data. 


5.2.2  Boundary  Conditions 


The  system  of  equations  is  to  be  solved  on  the  region 
® ^ walls  are  considered  impenetrable  by  solid 

particles,  so; 


u = 0 at  X = 0,  X 


(5.14) 


For  the  gas  motion,  two  options  are  available.  If  the  boundary 
is  considered  impermeable,  the  requirement  is: 

u^  = 0 at  the  boundary,  or,  using  (5.4)  and  (5.14), 


•^  = 0 at  the  boundary  (x  = 0 or  X) 


(5.15) 


If,  on  the  other  hand,  the  boundary  represents  a screen 
or  filter  through  which  the  gas  can  pass,  the  condition  is: 


P = P„ (t)  at  the  boundary 

D 


(5.16) 


where  P„(t)  is  a specified  boundary  gas  pressure  which  may 

D 

vary  with  time. 


5.2.3  Transformed  Relations 

It  is  useful  to  use  the  above  determinate  expression 
(Eqs.  (5. 5) - (5 . 13 ) ) to  rewrite  the  fundamental  field  equations 
in  a convenient  way.  The  solid  mass  relation  is  unchanged; 


II  (On)  . 0 


(5.17) 


For  the  gas  mass  conservation  constraint  (5.2),  we  first  use 
(5.4)  to  eliminate  explicit  dependence  on  u^: 


It  I(l-6)P£l  t [(l-9)p,ul  = P,  II 


(5.18) 
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If  the  gas  is  isothermal,  Eqs.  (5.7a)  and  (5.8a)  yield  for 
the  term  on  the  right; 

3 

In  the  adiabatic  case,  the  term  is  (using  5.7b  and  5.8b); 


k 

^ (Y-1)C.,T  p 


(5 


3x 


(y-i)c„t 
V 0 


(5 


Using  a little  algebra,  it  is  possible  to  show  that  this  be- 
comes ; 


It  * h ta-e)p,ul  = 


a 6 3 
2 ^ 


(5 


where 


k(0,r  ) (y-1)C,-T  r^ 

a E (y-1)C-,T  S—  - ^ ° ° 


n 


X k(6) 


(5 


fO 


'fO 


and  where 

6=0  for  the  isothermal  case 
= 1 for  the  adiabatic  case. 

The  solid  momentum  equation  may  now  be  written; 


(5 


3 (Gu)  ^ 
~5t  3x 


(euu) 


where  K = (y-DC^T^/P^ 


eg  - 

3 

3x 

^ 'v 

3u 

3x 

const. 

Of) 


6 apf 

3x 


(5 

(5 


.19a) 

.19b) 

.20) 

.21) 

.22) 


.23) 

.24) 
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and  where,  as  before. 


6=0  for  the  isothermal  case 
= 1 for  the  adiabatic  case. 

The  boundary  conditions  are,  on  solid  velocity  (Eq.  (5.23)) 


u=0atx=0. 


(5.25) 


At  boundaries  which  are  impermeable  to  gas,  the  condition  on 
Eq.  (5.20)  is  just; 


(Pf)  = 0 at  the  boundary 


(5.26) 


At  "prescribed  gas  pressure"  (permeable)  boundaries,  the  con- 


dition is 


PR(t) 

o2  = 2 

Of 


HM) 


(Y-l) 


(5.27) 


5.2.4  Summary  of  Principle  Equations 

To  repeat,  the  entire  system  is  now  reduced  to  the 
following ; 

1.  Solid  mass  conservation: 


30  9 I a \ r\ 

+ 5^  (eu)  = 0 


(5.28) 


Gas  mass  conservation: 


It  t(i-6)P£l  t [(i-e)p,ui 


6 3 ot  9 / 2 \ 

~ 9x  2 3x  (^f) 


(5.29) 


Momentum  conservation: 


3t 


0 

u 

Pf 

g 

a 

Y 

6f 


K 

V 

6 


^fO 

r 


f y 
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(0u) 


+ 


3x 


3x 


= 1 - porosity 

= solid  velocity 

= gas  interstitial  density 

= gravity  acceleration 

- >:<6)  'V-DCyT/ynjo 

= ratio  of  specific  heats 

= ^ (6-0  ) ^ for  6 > 0 
* 0 0 

= 0 for  0 < 0 
— 0 

= y-  (x  + I p ) X (6/6  ) 

Pg  \ 0 3 0/  0 

= 0 for  isothermal  gas 

= 1 for  adiabatic  gas 

=:  gas  constant  volume  specific  heat 

= standard  gas  temperature 

= standard  gas  density  (at  T ) 

0 

= standard  gas  viscosity  (at  T ) 

0 

= solid  particle  radius 

= solid  particle  density 

= solid  "viscosity  coefficients"  at  0 


exp 


k(6) 


2(1-6) 
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(-4.093  + O.3625[l-0]) 


Thus,  to  specify  a problem,  the  following  must  be  provided: 
Gas  Properties:  6,  y,  C^,  T^,  Hjq/  P^q 

Solid  Properties:  / 6 / a^ , X , y , r 

So  0 0 0 

Domain:  g,  length  of  domain  (X) 


Boundary  Conditions:  P, (t)  and/or  P„(t)  (gas  pressures  at 

Li  K 

X = 0,  X = X) , or  specification  of  im- 
permeable boundary 

Initial  Conditions:  Distributions  of  6 (x) , u(x),  p^  (x)  at 

t = 0 


5.3  FLUB  CALCULATIONS  - SAND  COLUMN  EXPERIMENTS  (100  psi 

OVERPRESSURE) 

When  the  blast  wave  from  an  explosion  passes  over  a 
porous  permeable  land  surface,  compressed  air  will  be  driven 
into  the  soil  and  will  diffuse  to  greater  and  greater  depths 
as  time  goes  on.  Once  the  above-surface  blast  pressure 
declines,  this  air  will  begin  to  flow  back  upward,  out  of 
the  soil.  If  the  soil  is  loose  and  the  blast  pressure  suf- 
ficiently great,  this  "surface  breathing"  may  carry  consider- 
able amounts  of  solid  material  aloft  as  the  compressed  air 
rushes  out  of  the  pores.  This  process  has  been  termed  "re- 

r 8 1 

verse  percolation"  by  Zernow,  et  al.  ^ 

To  evaluate  this  effect,  Zernow,  ^ al . conducted  a 
series  of  small-scale  experiments  in  which  a vertical  shock 
tube  was  used  to  generate  a long-period,  100  psi  (6.8  atm) 
overpressure  air  shock.  This  shock  impinged  from  above  upon 
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a vertical  tube  partially  filled  with  sand;  the  diameter  of 
the  sand  column  was  8 inches  (20  cm)  and  its  height  was 
about  71  inches  (180  cm) . Measurements  were  made  of  the 
variation  with  time  of  the  position  of  the  sand  surface  and 
of  the  interstitial  gas  pressure  time-history  at  various 
depths  below  the  original  sand  surface  level.  These  experi- 
ments were  carried  out  with  various  shock  wave  periods  and 
with  various  "batches"  of  sand  characterized  by  different 
particle  sizes.  In  most  of  the  experiments,  the  bottom  of 
the  sand  column  was  "partially  vented";  small  holes  in  the 
bottom  of  the  tube  containing  the  sand  were  provided  to  allow 
restricted  air  passage.  A few  experiments  were  also  per- 
formed with  these  holes  blocked  off  ("unvented")  so  as  to 
completely  prohibit  air  motion  at  the  bottom  of  the  bed. 

The  observed  sand  lofting  was  much  greater  in  the  latter 
case. 

This  experimental  configuration  was  treated  using 
the  original  FLUB  one-dimensional  code,  both  to  test  the  code 
and  to  gain  insight  into  the  important  mechanisms  operative 
in  the  sand-lofting  process.  Two  complete  calculations  were 
made.  These  correspond  to  tests  using  the  "1.36  inch  dia- 
meter orifice"  and  "Batch  II  Standard  Sand"  (see  Zernow, 

r g I 

et  al. ) . '■  ^ The  experimental  air-pressure-versus-time  his- 
tory at  the  surface  was  fit  with  the  following  function: 

P(0,t)  = 1 + 6.8  f (t)  exp  (-t/0.21  sec) 

where  P(0,t)  is  the  absolute  pressure  at  the  surface  in 
atmospheres,  and 

f(t)  =1  (t  < 0.55  sec) 

= 1-4  (t-0.55)  (0.55  sec  f.  t £ 0.8  sec) 

=0  (t  > 0.8  sec) 
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This  corresponds  to  a surface  burst  of  1.0  MT.  The  sand  bed 
was  represented  by  72  Eulerian  zones,  each  of  height  2.5  cm 
for  a total  height  of  180  cm;  an  additional  28  empty  zones 
were  provided  above  the  bed  to  permit  up  to  70  cm  of  sand 
lofting  before  material  strikes  the  upper  mesh  boundary. 

A total  of  360  Lagrangian  nodes  were  used  to  follow  sand 
motion. 

The  air  was  treated  as  an  adiabatic,  y = 1.4  ideal  gas; 
the  pressure/volume  relation  was  the  isentrope  reached  by  am- 
bient air  shocked  to  6.8  atm.  The  increase  in  air  viscosity 
with  temperature  was  taken  into  account.  Most  of  the  essen- 
tial sand  properties  are  available  in  Zernow,  ^ al . ; the 

crystal  density  was  2.4  g/cm^ , the  initial  porosity  (<))  ) was 

0 

taken  as  uniform  throughout  the  bed  and  equal  to  0.33,  and 
the  particle  diameter  (required  by  the  code  for  permeability 
calculations)  was  taken  as  650  microns,  which  is  given  by 
Zernow,  et  al . as  the  mean  particle  size  for  "Batch  II  Standard 
Sand" . The  bulk  pressure  in  the  sand  was  taken  as 

Pgg  = 3 ^1  - J ^ X 10®  atm  for  0 < 4)^  = 0.33 
0 

= 0 for  ())  > (j) 

0 

(where  <}>  = porosity)  which  agrees  reasonably  well  with  the  ex- 
perimental measurements  over  the  pressure  range  of  interest. 
Zernow,  ^ al . did  not  measure  the  intrinsic  viscosity  of  the 
sand  bed;  for  the  calculation,  the  viscosity  was  taken  as 
zero.  Preliminary  calculations  showed  the  results  to  be  in- 
sensitive to  reasonable  values  of  the  parameter. 

Two  cases  were  run;  the  first  treated  the  bottom  of  the 
computing  region  as  completely  impermeable,  corresponding  to 
the  experimental  "unvented"  case.  The  second  allowed  complete 
air  mobility  at  the  bottom  of  the  mesh,  so  that  the  gas  pres- 
sure there  remained  at  one  atmosphere  throughout  the  calcula- 
tion. This  "fully  vented"  calculation  was  intended  to  simulate 
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the  "partially  vented"  experimental  case,  but  of  course  is 
not  quite  the  same  thing.  Exact  simulation  of  the  experi- 
mental boundary  conditions  would  have  involved  a fair  amount 
of  code  development,  and  was  therefore  not  attempted. 

The  height-versus-time  of  the  surfaces  of  the  sand 
columns  computed  by  the  FLUB  code  are  compared  with  measure- 
ments in  Figure  5.2.  The  early  "precursor"  cloud  in  the  ex- 
perimental measurements  is  associated  with  high  frequency 
noise  in  the  incident  pressure  signal  resulting  from  shock 
reverberations  in  the  apparatus  as  Zernow,  ^ al.  have  demon- 
strated experimentally.  Measured  peak  heights  for  the 
partially  vented  case  varied  between  16  and  19  cm;  the  com- 
puted result  for  the  fully  vented  case  was  15.93  cm.  For  the 
unvented  case,  the  experiments  were  unable  to  determine  peak 
height,  but  it  was  known  to  be  > 35  cm;  the  computed  peak 
height  was  46.58  cm.  The  agreement  is  seen  to  be  quite  good, 
particularly  in  view  of  the  experimental  test-to-test  varia- 
tions and  the  inability  of  the  code  to  exactly  simulate  the 
lower  boundary  condition  for  the  vented  case.  Interstitial 
gas-pressure  histories  at  depths  of  zero,  84  and  146  cm  are 
shown  in  Figure  5.3.  For  comparison,  the  experimental  results 
for  the  partially  vented  case  are  listed  below  (no  results 
were  presented  by  Zernow,  et  al . for  the  unvented  case) ; 

Depth 

84  cm  146  cm 

Peak  gas  pressure  (atm)  1.34-1.71  ''^0.6 

Time  of  peak  (sec)  0.22  - 0.23  ^<0.3 

Again,  agreement  is  good. 

The  total  mass  of  sand  above  the  pretest  surface  per 
unit  area  as  computed  by  the  code  is  shown  as  a function  of 
time  for  each  case  in  Figure  5.4.  As  can  be  seen,  the  un- 
vented case  lofts  substantially  more  material. 


i FLUB  computations 


Gas  Overpressu 


Time,  msec 


Figure  5.3.  Gas  overpressure  versus  time  at  surface  and 
z = 84  and  146  cm  as  computed  by  FLUB  code. 


Figure  5.4.  Total  solid  mass  aloft  versus  time  as  computed  by  FLUB  code 


5.3.1  Remarks 


Based  on  the  results  in  Figure  5.4,  it  would  be 
reasonable  to  assume  that  at  least  10  g/cm^  would  be  avail- 
able for  entrainment  in  the  negative  phase  af teirv/inds.  If 
we  assume  that  this  is  a representative  mass  lofting  poten- 
tial, for  overpressures  greater  than  100  psi,  then  for  a 
1 MT  surface  burst  (100  psi  peak  overpressure  occurs  at  a 
range  of  1 km)  approximately  0.3  MT  will  be  lofted.  This 
corresponds,  roughly,  to  the  predicted  amount  of  dust  that 
is  sheared  off  the  surface  in  a typical  1 MT  SCOUR  calcula- 
tion (see  Figure  3.9). 

It  is  important  to  stress  that  the  FLUB  predictions 
indicate  that  the  dust  is  lofted  during  the  transition  to 
the  negative  phase  and  that  the  dust  is  not  violently  "pop- 
corned"  into  the  free  stream.  Rather,  the  surface  dust  has 
been  fluidized  down  to  depths  of  approximately  200  cm. 
Typical  dust  concentrations  in  this  fluidized  layer  are 
0.5  g/cc.  Thus,  under  the  dry,  incohesive  soil  conditions 
evaluated  in  both  the  experiments  and  FLUB  calculations,  it 
would  appear  that  the  surface  boundary  conditions  could  be 
significantly  different  than  that  assumed  in  the  SCOUR  shear 
lofting  model.  Instead  of  an  undisturbed  dust  surface,  over 
which  a saltation  layer  develops  in  the  negative  phase,  the 
early  negative  phase  winds  would  entrain  dust  from  a fluid- 
ized layer  as  demonstrated  in  Figure  5.5. 


Equally  important  to  remember  is  that  the  SCOUR  model 
for  dust  lofting  is  quantified  through  a mass  lofting  rela- 
tion 


ND 


V 


<<>a“5  - 


(5.31) 


where  M^  is  the  aerodynamic  lofting  rate  of  dust  from  the 
particulate  surface  and  the  shear  (p^u^)  is  determined  by  a 
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(Air  Velocity) 


Figure  5.5.  Schematic  of  new  air-soil  interface  that  results  if  surface  soil  is 
fluidized  during  reverse  percolation. 


rough-wall  shear  law  for  the  so-called  "perfectly  rough" 
flow  regime, 

u^  ^ 0.707 

^e  ^2.87  + 0.7  Hn 


(5.32) 


Hence,  if  one  may  continue  to  assume  that  the  boundary  layer 
is  still  in  the  perfectly  rough  regime,  and  that  an  effective 
(rough)  surface  is  buried  in  the  fluidized  layer,  then  the 
entrainment  at  the  surface  may  still  be  proportional  to  the 
shear  at  the  fluidized  layer-air  interface  (as  postulated  in 
Eq.  (5.31)).  The  SCOUR  model,  if  all  these  hypotheses  are 
appropriate,  would  still  provide  a reasonable  means  of  cal- 
culating the  negative  phase  mass  lofting  potential  of  nuclear 
blast  wave  flow  fields.  Presently,  a minimum  roughness  height 
of  10  cm  has  been  incorporated  in  SCOUR  to  take  this  fluidiza- 
tion effect  (and  surface  modulations)  into  account. 

5.4  FLUB  CALCULATION  - 1100  PSI  OVERPRESSURE 

The  surprising  degree  of  surface  soil  fluidization 

predicted  at  the  100  psi  overpressure  provided  the  impetus 

for  additional  FLUB  calculations  at  higher  overpressures. 

A representative  case  is  the  1100  p.si  overpressure  station 

for  a 36.6  KT  detonation  at  a 350  ft  HOB.  This  corresponds 

to  typical  yields  and  HOB  for  NTS  events.  To  better  match 

NTS  conditions,  particle  size  was  reduced  to  0.1  mm  particle 

diameter  and  soil  porosity  to  0.25.  To  maximize  the  reverse 

percolation  effect,  the  "unvented"  boundary  condition  (at  a 

depth  of  180  cm  below  the  ground  surface)  was  assumed.  The 

overpressure  time  history,  graphically  depicted  in  Figure 

[37] 

5.6,  was  computed  from  Brode's  analytic  expression  for 

a range  of  165  meters.  As  indicated,  the  positive  phase  is 
complete  0.3  seconds  after  air  shock  arrival.  All  other 
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Figure  5.6.  Overp 


aspects  of  the  calculation  were  identical  to  the  pareuneters 
and  grid  used  in  the  preceding  100  psi  test  case. 


The  results  for  this  FLUB  calculation  are  presented 
in  Figure  5.7  (soil  altitude),  Figure  5.8  (blowoff  velocity), 
and  Figure  5.9  (mass  lofted).  As  expected  on  the  basis  of 
the  100  psi  results,  the  unvented  soil  column  was  almost 
completely  fluidized  by  the  end  of  the  positive  phase,  and 
the  calculation  could  not  be  continued  into  the  negative 
phase  without  a major  effort  to  rezone  the  problem.  The 
conclusions  cited  in  Section  5.3.3  apply  to  these  results; 
i.e.,  reverse  percolation  at  high  overpressures  in  cohesion- 
less, porous  soils,  is  actually  a fluidization  process  that 
may  be  responsible  for  the  establishment  of  a new  boundary 
condition  for  surface  dust  entrainment  in  the  early  negative 
phase.  A major  question  not  answered  in  this  calculation, 
due  to  grid  size  limitations,  was  the  time  for  the  fluidiza- 
tion process  to  decay.  Secondly,  since  blowoff  velocities 
are,  at  most,  the  same  order  of  magnitude  as  the  near-surface 
air  velocities,  the  interaction  between  the  subsurface  (per- 
colation) air  flow  and  the  wind  field  will  seriously  affect 
the  development  of  the  fluidized  bed. 

It  should  be  emphasized  once  again  that  the  present 
results  assume  that  the  soil  has  no  cohesion  (relative  to  the 
pressure  gradients  imposed  by  the  air  in  the  soil  pores) ; that 
the  porosity  is  relatively  high  (25  percent)  and  is  distri- 
buted uniformly  throughout  the  soil  surface;  and  that  there 
are  no  shock  waves  in  the  soil  matrix.  Moreover,  because 
of  the  surface  wind  interaction  effects,  the  problem  should 
be  considered  in  two  space  dimension.  Nevertheless,  these 
preliminary  results  give  evidence  that  the  general  problem  of 
(air)  boundary  layer  interaction  with  a fluidized  surface 
layer  should  be  more  closely  studied. 


112 


Figure  5.8.  Upward  velocity  of  topmost  particles  in  fluidized  surface  region,  from 

FLUB  calculations  for  1100  psi  overpressure  station,  36.6  KT  at  300  ft  HOB 


Figure  5.9.  Mass  of  dust  in  fluidized  surface  layer  above  the  original  surface,  from 

FLUB  calculations  for  1100  psi  overpressure  station,  36.6  KT  at  300  ft  HOB 
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APPENDIX  A 
SC0UR2  INPUT 


Card  input  for  SC0UR2  contains  the  basic  parameters 
necessary  for  a ECOUR  calculation.  All  units  are  cgs. 


1. 


RUN  Parameters 

CYCLE 

= 

cycle  to  begin  on. 

TDUMP 

= 

time  to  begin  on. 

YIELD 

= 

yield  in  megatons. 

HOB 

= 

height  of  burst. 

WC 

= 

soil  water  content. 

mass  fraction. 

TND 

= 

time  to  end  on. 

EDIT  Parameters 

NPRINT 

= 

scale  factor  cycles 

between 

zone  prints. 

NDUMP 

= 

scale  factor  cycles 
dumps. 

between 

restart  tape 

NPPRT 

= 

scale  factor  cycles 
prints. 

between 

particle 

NCPRT 

= 

scale  factor  cycles 

between 

printer 

plots , 


GRID  Parameters 


iriAX 

JMAX 

DXZ 

DYZ 

KEMAX 

PEMAX 

TEMAX 


= number  of  horizontal  zones. 

= number  of  vertical  zones. 

= zone  width  (constant  in  radial  direction) , 

= scaling  factor  for  concentration  (printer 
plots) . 

= scale  factor  for  kinetic  energy  printer 
plots . 

= scale  factor  for  potential  energy  printer 
plots . 

= scale  factor  for  total  energy  printer 
plots . 


121 


DPriAX 

CCMAX 

DPCMAX 

RHSOIL 

RHOP 

GRAY 

GAMMA 

XSTART 


= scale  factor  for  dynamic  pressure 
printer  plots. 

= scale  factor  for  concentration  ratio 
printer  plots. 

= scale  factor  for  dynamic  pressure  ratio 
printer  plots. 

= density  of  soil. 

= crystal  density  of  soil  particles. 

= acceleration  due  to  gravity  (positive  is 
downward) . 

= time  integration  coefficient,  shock 
lofting  integral  at  shock  front. 

= range  of  first  zone  beyond  crater. 


No.  of  particles  = number  of  particle  sizes  in  the 

calculation. 

Initial  random  number  = initial  number  for  random 

number  subroutine. 


PARTICLE  Parameters 
RCHAR 


RCHARN 
RADZ (ND) 
IFRONT 


characteristic  particle  radius  for  which 
the  roughness  height  (HBAR(I))  is  ca] 
culated  (positive  phase) . 

RCIiAR  (negative  phase)  . 

particle  radius  (for  each  size  particle) . 

flag  for  creating  first  particle  of  size 
ND  at  shock  front  (0  = no,  1 = yes) . 


fiASS  LOFTING  Parameters 

DEAR  = minimum  roughness  height. 

FACN  = boundary  layer  diffusion  parameter  for 

the  negative  phase. 

FACP  = boundary  layer  diffusion  parameter  for 

the  positive  phase. 

SMALLA  = fraction  of  mass  of  shock  lofted  SCOUR 

particle  necessary  to  create  wind-lofted 
SCOUR  particle. 

DELMAX  = maximum  altitude  of  the  boundary  layer. 


ALPHA (ND)  = scale  factor  in  prescription  of  particle 
ejection  velocity,  particle  size  ND. 

SMALLB(ND)  = mass  fraction  of  particle  size  ND. 

XLAJl(ND)  = mass  lofting  coefficient  of  particle 

size  ND. 

7.  CRATER  LOFTING  Parameters 

PLM  = mass  represented  by  SCOUR  crater  particles. 

RCR  = crater  radius. 

DCR  = crater  depth. 

FORM  = fraction  of  crater  mass  to  be  ejected. 

YLOAD  = maximum  ejection  altitude  of  crater 

particles  at  TDUMP. 

VBAR  = crater  particle  ejection  velocity. 

8.  GLAB  = alphanumeric  label  for  the  run. 


SC0UR2  OUTPUT 


Printout  of  the  Input  Parameters  (Typical) 


RUN  PARAMETERS 

YIELD*  l•000*00 
HOB*  o<aoo 
WC*  0.000 
TNO*  I.e00*02 

EDIT  PAKAME1ER5 
NPHINT*  10 
NDUMP>  10 
NPPRT*  10 
NCPRT*  10 

GRID  parameters 


IMAX* 

75 

JMAX* 

20 

DX2* 

5.C00»03 

DT2* 

1 .000-03 

KEMAX* 

1 .000-03 

PEMAX • 

1 .000-03 

TEMAX* 

1 .000-03 

OPMAX* 

1 .000*05 

CCMAX* 

1 .000-03 

OPCHAX* 

1 .000*05 

rmsoil* 

1 . 800*00 

RHOP* 

2.500.00 

GRAV* 

9.800*02 

GAMMA* 

1 .000*00 

xstart* 

1 . 750*09 

ANGLE* 

1 .000*00 

NO.  OF  PARTICLES*  I 

initial  RANDOM  NUMBCP* 

PARTICLE  parameters 
RCHAR*  3.S0O-O3 
RCHARN*  R. 000*02 
RAD2*  S.OCO-03 
IFRONT*  I 

MASS  L0FTIN6  PARAMETERS 


OBARa 

1.000*01 

FACf4« 

8.000-02 

F ACP* 

8.000-02 

SMALL  A* 

2.C00-01 

D£LHAX« 

1 .750*09 

ALPHA- 

1 .000*00 

SHALLB* 

9.000-01 

XL  AM* 

1 .000-01 
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APPLE  27 

c«ater  lofting  PAHAHETERS 

PLH-  i.ooo*ce 

RCRa  I.TSO^OR 

OCR*  8.7SO*C3  > 

FCRH*  2.000-02  ! 

TLOAD*  I. 000^02  I 

VBAR*  1.000*08 

crater  LOFTInG:  crater  HASS-  ».8708*ll  NO.  OF  PRELOAOEO  PARTICLES-  |89 

Crater  mass  is  the  mass  of  dust  that  occupied  the 
crater.  The  number  of  preloaded  particles  refers  to  the 
number  of  SC0UR2  particles  used  to  represent  the  fraction 
of  the  crater  mass  that  is  to  be  tracked  in  the  SCOUR2  cal- 
culation. The  APPLE27  is  the  alphanumeric  label  for  this 
particular  run.  Note  that  in  this  calculation  only  one 

particle  size  (radius  = 5 x 10  ^ cm)  was  assumed  for  the 

soil,  and  that  since  RCHAR,  RCHARN  were  specified  as  dif- 
ferent values,  the  characteristic  roughness  height  (h)  was 
not  dependent  on  the  particle  altitudes  (h  = DEAR) . Further- 
more, the  radius  of  the  crater  particles  is  missing  from  this 
printout.  It  will  be  included  in  subsequent  SC0UR2  input 
printouts . 
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2. 

APPLC  21 


BESIAVAIUBIE  COPY 


Boundary  Layer  Information 

!•  l.Oltl'OO  IS*  I.Sflf^oS 


1 

U2 

OCLTA 

HOAR 

voirr 

ustar 

RHOZ 

V2 

eotPTM 

5UM  fOTSuM 

2 

1 

•S0OI*O9 

4.6370*02 

1.0000*01 

4 

«4|3t*00 

2 

•0569*03 

3.9290-04 

1 .4353*04 

2.3744-02 

2.6854*07 

J 

1 

*SS0I*09 

4.2874*02 

I .0000*01 

4 

•9409*00 

1 

•4937*03 

3.7784-04 

•1 .3945*04 

1 .8954-02 

?.4C69*07 

H 

-1 

• 09B2*0*4 

J.Bl 13*02 

1.0000*01 

3 

• 1 827*00 

1 

•2321*03 

3.B7V2-06 

•1 .5358*04 

2.9301-02 

9.9708*07 

S 

• 1 

• 1 1 7*4*o'4 

3.2287*02 

1 .0300*01 

2 

•Bo72*00 

1 

•2966*03 

4.1291 -04 

4.5200*03 

3.7545-02 

7.9309*07 

4 

-2 

• 1728*03 

9.5257*01 

1.0030*01 

1 

• 4490'0 1 

3 

.1041 *02 

4.2452-06 

5.81 1 1*03 

2.2409-02 

5. 1 582*07 

7 

2 

•6943*03 

t . 1 1 30*02 

1 •0500*01 

2 

.4495-01 

3 

.740  1 *C2 

4 .390 1 -04 

4.4268*03 

2 .9969-02 

7.6243*07 

• 

1 

•S7d3*03 

7 .5909*01 

1 .0000*01 

1 

•0024-01 

2 

•3430*02 

4 .6682-04 

3.2414*03 

4 .705  1 -(J2 

1.3)03*08 

9 

2 

• 304«4  *02 

2.1701*01 

1 *0000*0  1 

4 

•2956-03 

4 

•8330*01 

5.6394-06 

2.6465*03 

5 . 22  10-02 

1 .62)8*08 

io 

8 

.‘4230*Oi 

1.0969*01 

1 «0C30*UI 

8 

• 1 395-04 

2 

•0753*01 

7.4431-06 

1.6165*03 

6 • 9656-02 

2.  |937*0« 

1 1 

-1 

.•4  7 7 9*03 

7.2965*01 

1 *0030*01 

9 

•4674-02 

2 

.2339*02 

9.9741 -04 

6.5896*02 

9.6529-02 

1 .7 10)«08 

12 

-1 

. 4 1 .>  1 *03 

9.0437*01 

1 •0000*ul 

1 

•3564-01 

2 

.3054*02 

1 .5332-05 

3.650 1 *02 

2 .90  12-02 

1.1601*08 

1 3 

- 1 

.9389*03 

8.7809*01 

1 .3000*0  1 

1 

.629C-01 

2 

.8170*02 

2.2000-05 

3 . 65 1 7*02 

7 • 5 746-02 

3.2125*08 

I** 

-8 

.•472S*02 

4.9995*01 

1 .0000*01 

4 

• 1 268-02 

1 

.3974*02 

2.8650-05 

3.C696*02 

5.7959-02 

2 .599  2 • 08 

IS 

-3 

• •4990*02 

4 • 5 1 6 1 *03 

1 .0300*01 

2 

.2767*00 

2 

. 5 1 98  *0 i 

3.6469-05 

2 .27  J 9*02 

6.9  1 29-0? 

3.  )225»C8 

14 

-1 

. 9939*02 

1 . 7500*04 

1 .0000*01 

3 

.4524*00 

1 

•4618*01 

5.0399-05 

1 .5329*02 

V.  1 942-02 

9,6559*08 

1 7 

-8 

.2*499  *02 

1 .0579*04 

1 .0000*01 

9 

.4336*00 

4 

•0437*01 

6 .707  3-05 

1 . 1519*02 

9.5o9V-02 

5. 1 992*08 

1 s 

• 6 

« 1 734*02 

1 .7500*34 

1 .0000*01 

1 

•2049*01 

4 

.5436*01 

6.5558-05 

1 .31  10*02 

5.6080-02 

3.1712*08 

1 » 

7 

• 1 390*02 

1 . 7500*04 

1 •000C*(3l 

1 

.444 1 *01 

5 

.27jC*0l 

1 .0043-04 

1 .9565*02 

8.6928-02 

5.1317*08 

20 

; 

. 3S29*02 

1 .7500*04 

1 .0000*01 

1 

•5406*01 

5 

.9450*01 

1 .9015-04 

8. 7Vti3*0l 

7.0939-02 

9,3612*08 

21 

A 

• I4oa7*02 

1 .7500*04 

1 .0000*01 

1 

.9272*01 

4 

.5037*ul 

1 . 8 1 J 3-OH 

6.9895*01 

9.2593-02 

4.0219.08 

22 

1 

•3211 *03 

1 .7500*04 

1 .0000*01 

2 

.9825*01 

9 

.B3b5*0t 

2.3 1 9H-04 

8 . 1 955*01 

8.5502-0? 

5.faC2l *38 

2 J 

1 

.41 33*03 

1 .7500*04 

1 .0000*01 

3 

.7691 *01 

1 

.?0‘««*02 

2.9251-04 

7.6396*01 

7.6601-02 

5.9  J 96*08 

2‘< 

I 

.9308*03 

1 .7500*04 

1 .0000*01 

4 

.7248*01 

1 

.4462*02 

3.6005-04 

6.8,^60*01 

6.5997-02 

9.6960*08 

2'i 

2 

. 2*400*03 

1.7500*04 

1 .0000*01 

5 

•7214*01 

1 

.4332*02 

4 . 3235-04 

6.9783*01 

1 . 0C«  7-0 1 

7. 7003*08 

2 

.S820*03 

1 .350 1 *04 

i *0000*0  1 

5 

.3209*01 

1 

•9467*02 

5.07  10-0't 

6.1671*01 

6.179  2-02, 

9 . 69  2U  • 04 

77 

2 

.9392*03 

9 • 9608*03 

I •0000*01 

9 

•6632*01 

2 

•2238*02 

5. 8232-09 

9.9902*01 

3.8541-02- 

3.  1 635  *08 

28 

3 

.31 IS*03 

7.4913*03 

1 •0003*01 

4 

. 1 4Q4*0I 

2 

.5147*02 

6.5420-09 

5.7736*01 

9 . 16  8 1-02 

7.7767*0 

29 

3 

• 7 1 10*03 

5.6844*03 

1 •0000*0  1 

3 

•7328*01 

2 

.8291*02 

7.2732-09 

5.7289*01 

5. 7066-02 

5.C02 • *08 

30 

•4 

• ) 5*44*03 

9.31  19*03 

1 .0000*01 

3 

.3476*01 

3 

• 1 803*02 

7.957  J-04 

5.7396*01 

6.6359-02 

6.00)8*0'’ 

31 

i 

.4037*03 

3.3276*03 

1 *3000*01 

3 

•0329*01 

3 

.5396*02 

a. 4023-09 

5.7973*01 

9.5996-02 

9 . 2906  *fl'1 

32 

s 

•0763*03 

2.5623*03 

1.0000*01 

2 

• 76CC*0 I 

3 

.9228*02 

9.2127-09 

6 . M 23*01 

9 . 6 350-92 

9 . 2 6 2 9 • C-  f 

33 

s 

• 6 144*03 

1 .9645*03 

1.0000*01 

2 

•504p*0l 

4 

.3404*02 

9.801 3-09 

6.2515*01 

5.7992-02 

5.7369.0" 

3H 

6 

. 2299*03 

1 • 5o78*o3 

1 «0CC0*01 

2 

•2426*01 

4 

. 8622*02 

1 .03  7 7-03 

6.5813*01 

6 . 3 3 36-02 

6.9H69.3A 

3S 

6 

.4779*03 

1 • 1479*03 

1 .0003*01 

2 

•049 1 *01 

5 

•3996*02 

1 .0928-03 

7.1179*01 

2.6729-02 

2. 7954*34 

34 

7 

.5599*03 

8.7303*02 

1 .0000*01 

•8585*01 

5 

.9727*02 

1 • 1 449*03 

7.3373*01 

4 . 1525-02 

9.9615*04 

37 

a 

.277S*C3 

4.6214*02 

1 .0000*01 

•4857*01 

4 

•5851 *02 

1 .2007-03 

7.9900*01 

5. 8299-02 

6.9287*08 

3B 

9 

•0327*03 

4.9895*02 

1.0030*01 

.5270*01 

7 

•2411 *02 

1 . 2544*03 

6.7210*01 

5.0108-02 

5.66  71  *08 

3» 

9 

•8329*03 

3.7  1 7 8*02 

1 .0000*01 

•3787*01 

7 

.9498*02 

1 . 3099-03 

9.0282*01 

5. 7950-02 

6. 7 1 79 

*»0 

.0461*0*4 

2.7257*02 

1 .0000*01 

•2381 *01 

B 

.7199*02 

1 . 3456-03 

9.3220*01 

6.5963-02 

7.4332*04 

9 1 

. 1 56 ( *0H 

1 .96 13*02 

1 .0000*01 

•1054*01 

9 

•5436*02 

1 .4237-03 

1 .0308*02 

3.9936-02 

9.2975*08 

92 

.2*4  1 9*0*4 

1 .3657*02 

1 .3000*01 

9 

•8246*00 

1 

•0306*03 

1 .461 3-03 

1 .0695*02 

2.6336-07 

3 . 2 7 6 9 * f,  8 

•43 

.3296*0*4 

9,5145*01 

1 • ocoo  * L 1 

e 

•4852*00 

1 

• 1 255*03 

1 .5382-03 

1.1135*02 

9 . 2250-02 

5. 3 756  *08 

HS 

.*420  7 *04 

4.1139*01 

1 .3000*01 

7 

.5443*00 

1 

.2309*03 

1 .4CJ10-03 

1.1914*02 

3 . JU  I 7-02 

9.2992*06 

•4S 

.5609*04 

3.4341*01 

I .0000*01 

4 

•2962*C0 

1 

.3862*03 

1 .4850-03 

1 .2920*02 

6.9b07-03 

9.2)68*07 

•44 

• 7992*04 

1 .5693*0  1 

1 .0000*01 

4 

•9446*00 

1 

•4130*03 

1 . 7999-0) 

1.3964*02 

J . 2 7 1 8-02 

1*  7 26  1 «Q8 

•47 

.6533*04 

5.0930*00 

1 .0030*01 

3 

.7553*00 

1 

.8967*03 

1 .8765-03 

1 .3807*02 

1 . 2975-02 

1 . 7976*08 

S8 

0 

• ecco 

n.tooo 

1 .0000*01 

0 

• OOOw 

0 

• OOOG 

0.0000 

0.0000 

0.03C0 

O.COOO 

•49 

c 

• ccoo 

b.GC'OO 

1 .OCCO^Ol 

0 

• OCOO 

0 

.0000 

0.0000 

0 .0000 

O.OOCO 

O.COOO 

SO 

0 

• CPOO 

c • ocoo 

1 .0000*01 

0 

• 0000 

0 

.0000 

0.0000 

o.cooo 

0 *0300 

O.COOO 

S 

0 

• crco 

O.DOOO 

1 .0000*01 

0 

• cooo 

0 

• OCOO 

0.0000 

0.0000 

0*0300 

O.OCOO 

T = time  (sec)  . 

XS  = shock  front  (cm) . 

UZ(I)  - horizontal  velocity  of  air  at  center  of  zone  i (cm/sec). 

DELTA(I)  - boundary  layer  height  for  zone  i (cm). 

HBAR(I)  - average  height  of  particles  in  zone  i (cm). 

VDIFF(I)  - extra  diffusional  velocity  due  to  boundary  layer 
growth  in  zone  i (cm/sec) . 

USTAR(I)  - shear  velocity,  zone  i (cm/sec). 

RHOZ(I)  - air  density  of  zone  i (g/cm^). 

VZ(I)  - vertical  velocity  of  air  at  center  of  zone  i (cm/sec). 

EDEPTH(I)  - erosion  depth,  zone  i (cm). 

SUMTOTSUM(I)  - total  mass  (g)  of  all  particle  sizes  lofted  from 
zone  i up  to  time  T. 
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AVAILABLE  COPY 


3. 

APPLE  2? 


Mass  Lofting  Information 


)• 

0393*00  XS« 

R2 

2*9919*0S 

ENCRZ 

tau 

VZCHOZ 

OOTnZ 

ZMASS 

SUM2 

TOTSUH 

2 

2 

C000*09 

|.S34S*12 

1 .4433*01 

1*9473*03 

1 *9972*05 

1 

*8434*07 

•2*4369*04 

2*4854*07 

3 

2 

S000*09 

I«S|C2*I2 

1 •OB9o«ol 

1 *2489*03 

2*9341*05 

1 

*8159*07 

•1*0512*07 

2*4069*07 

3 

0000*09 

1 *9992*12 

S*8886*00 

1 * 1228*03 

1 *4919*05 

3 

*9241*07 

• 1 • 1844*07 

9*9708*07 

& 

3 

SC00*09 

1 *9984*12 

4*9351*00 

9*0290*02 

2*891 1 *05 

9 

*9718*07 

«| *9934*07 

7.9309*07 

4 

9 

C000*09 

1 *9794*  12 

9*1  iSo-Ol 

2*9936*02 

0*0030 

3 

•9354*07 

2*2402*04 

5*  1582*07 

7 

9 

S00C*09 

1 *9  150*1 2 

4*191C*0I 

1*9308*02 

0*0000 

9 

*9278*07 

7*9302*04 

7*4243*07 

8 

S 

CQC0*09 

1 *2909*1 2 

2*7294-01 

2*1658*C2 

0*0000 

7 

*2931 *07 

•2*3977*07 

1 .3303*08 

f 

s 

S00C*09 

1 .0772*12 

1*3623*02 

9*2434*01 

0*0000 

* 1 8 76*08 

•4.9 1 30*07 

1 *4238*08 

10 

4 

C0JO*O9 

8. 1692*1  1 

3*2055-03 

1 *0885*01 

Q.COOO 

*2409*08 

•9.2209*07 

2*1937*08 

1 1 

4 

SC00*09 

S. 7173*1  1 

9*9773-01 

1 *6928*02 

o.ooco 

7 

*9346*07 

1 *0996*04 

1 *7103*08 

12 

7 

0000*09 

3*8995*1 1 

8. 1989-01 

2*0099*02 

0*0000 

6 

*35U2*07 

-1 .2729*08 

1*1601 *C8 

13 

7 

S000*09 

2*4703* 1 1 

1 * 7959*00 

1*5995*02 

9*9088*05 

*3519*08 

•1 *9963*07 

3*2125*04 

M 

8 

0000*09 

1 .8124*1 1 

5*5996-01 

1*2032*02 

0*0000 

*2589*08 

•1 *2992*08 

2*5992*04 

IS 

B 

So00*09 

1 * 1949*1  1 

2*3156-02 

1*9393*01 

O.OCOO 

*9927*00 

*•2  *9648*07 

3.  3225*0H 

14 

9 

C000*09 

7.91  15*10 

1 *0749-02 

1 *2613*01 

o.ocoo 

•7093*00 

-2.01C9*07 

9*6559*08 

17 

9 

S000*09 

S .998  1 * 1 0 

2*9956-01 

9.5899*01 

0*0000 

2 

* 2065*08 

•2 • 786H«0Q 

5.  1992*08 

18 

oooo*os 

3.7849*10 

1 .7702-01 

2*4667*01 

o.ocoo 

2 

•9052*08 

-2* 1008*06 

3*1712*08 

19 

csoc*os 

2*40?0*I0 

3.0150-01 

3*9606*01 

o.ococ 

* 98 1 9*06 

• 9.  7580*07 

S. 1 31 7*08 

20 

iooo*os 

1 .8221  * 10 

9*1553-01 

3*3800*01 

0*0000 

*9233*06 

• 1 .0408*08 

9 • 3a 1 2 *04 

21 

IS30*0S 

1 .31 18*10 

7*851 l-Ol 

5*9586*01 

c.coco 

*0916*03 

•2  . 1 058*07 

4*0219*38 

22 

2noo*os 

9.7670*09 

2.2996*00 

8*B|98*01 

2*4598*04 

2 

* 1 929*00 

-9.6254*00 

5 • 8021*08 

23 

2500*OS 

7.5607*09 

9 * 2959*00 

7.935c*01 

8*5721 *04 

*9037*08 

-7.5219*07 

5*  9 1 96*08 

2‘* 

3000*0s 

4.0459*09 

7*5308*00 

8*955| *01 

1 .9892*07 

•5992*08 

•1 . 1915*09 

9 . 8980*08 

2S 

3500*0S 

5.0325*09 

1.2299*01 

1 *0741  *02 

2*2144*07 

9 

. 1376*08 

-1 .4908*08 

7*7003*08 

24 

93U0*0S 

9*3037*09 

1 *9218*01 

1*9518*02 

2.7595*07 

*51 76*00 

•6.9945*07 

9.6920*08 

27 

9SOO*OS 

3*7812*09 

2*6796*01 

1*17  16*02 

5*9035*07 

.6739*00 

-1 .01 13*00 

3.  1 635*08 

28 

sooo*os 

3.9307*09 

9 * 1996*0 1 

2.3197*02 

9*1  139*07 

2 

•0768*08 

•2. 1 309*08 

7.7767*08 

29 

ssor*os 

3.  12S9*09 

5*8213*01 

2.9295*02 

5*7955*07 

7 

* 1 335*07 

•2.0599*07 

5*002 1 *08 

30 

4030*0S 

2.9206*09 

6*0989*01 

2*7939*02 

7*2802*07 

•9990*08 

•5.3335«07 

6.003&*C8 

31 

6SOO*OS 

2.7473*09 

1 .0770*02 

3*2135*02 

8.4121*07 

2 

. 1 729*00 

-9.4962*07 

9*2906*08 

32 

7aoo*os 

2.4516*09 

1*9177*02 

3*1203*02 

1 *2097*08 

*5043*08 

9 . 1596*04 

9*2629*08 

33 

7530*0S 

2*5676*09 

1 *8635*02 

3*  1 728*02 

1 *6059*08 

•5009*08 

O.OCOO 

5*7389*08 

39 

B000*0S 

2.5079*09 

2*9532*02 

3*5698*02 

1*9352*08 

•9228*00 

0.0000 

6 • 9 9 69*08 

3S 

8SOO*OS 

2*9647*09 

3.1841*02 

5*2939*02 

1 *7603*C8 

• 1 962*08 

•2*963 1 *08 

2*7956*08 

36 

9CCO*OS 

2.9394*09 

9*0915*02 

3*5194*02 

3*9410*08 

2 

•1289*08 

• 7 . 8 1 1 1 * 0 > 

9*9615*08 

37 

9SOO*OS 

2*9236*09 

5*2044*02 

9*9928*02 

3*1601 *06 

•5592*06 

•1*1 399*00 

6*9287*08 

34 

2 

oooo*os 

2*9|68*o9 

4*5789*02 

9*7222*02 

9*3699*08 

•2924*08 

-5*7343*07 

5*667 1 *00 

39 

2 

DS00*OS 

2*9  1 78*09 

8*2759*02 

5*3594*02 

9*9659*00 

2 

.2701*08 

o.cooo 

6.7  1 79*08 

90 

2 

1 0J0*0S 

2.9256*09 

I *0385*03 

5*1551*02 

4*4308*08 

• 1 479*08 

•2*3738*08 

7.6332*08 

9 

2 

lS0O*0S 

2.9389*09 

1*2967*03 

5*7232*02 

7*6957*08 

.5137*06 

•2*7338*08 

9.2975*08 

92 

2 

2O00*0S 

2.9553*09 

1*5970*03 

1 *0215*03 

5*9019*08 

*0751 *08 

•2*4975*07 

3.2769*08 

93 

2 

2sao*os 

2.9730*09 

1 *9989*03 

4*2964*02 

1 *0931  *09 

*9673*08 

-3*9089*08 

5.3756*08 

99 

2 

3000*0S 

2.9982*09 

2*9257*03 

1*2259*03 

7*  1 956*08 

* 7392*00 

•2. 5551 *08 

9 . 2992*08 

MS 

2 

3SJ0*0S 

2.5344*09 

3*2973*03 

1 *2370*03 

7*7998*08 

9 

*23o8*07 

0*0000 

9.2368*07 

94 

2 

90C0*0S 

2*5888*09 

9*6685*03 

9*7200*02 

1 *9902*09 

•7241 *08 

0.0000 

1 .726  1 *06 

97 

2 

9S00*0S 

2*6222*09 

6.3995*03 

1 *304 1*03 

1 *5C03*09 

*7974*08 

c.oooo 

1 . 7976*08 

98 

2 

sooo*os 

0.30D0 

o.ocoo 

o.cooc 

0*0C00 

0 

*0000 

o.cooo 

0.0000 

99 

2 

S5J0*0S 

o.ooco 

o.occo 

0*0000 

0*0300 

0 

• 0000 

0.0000 

o.ocoo 

SO 

2 

4030*08 

0.0000 

0*0000 

o.occo 

0*0030 

0 

• 0000 

o.cooo 

o.ocoo 

SI 

T 

2 

6S00*0S 

tine 

0.0000 

(sec)  . 

O.OCOO 

0*0000 

0*0000 

0 

• 0000 

o.cooo 

0.0000 

XS  = shock 
R2(I) 

ENERZ (I) 

TAU ( I ) 

VZEROZ (I,ND) * 


front 


location  (cm) . 

X coordinate  of  center  of  zone  i (cm) 
energy  of  air  in  zone  i (ergs/^) 


ND 


ZMASS (I,ND) * 
TOTSUM(I,ND) * 


zone  1 ai 
particle 


size 


tide 

size 

ND 

in 

zone 

ND 

in 

zone 

boundary  layer  shear  (dynes/cm^). 
initial  velocity  of  lofted  particle  size 
in  zone  i (cm/sec) . 

DOTMZ(I,ND)*  - mass  lofting  rate  for 
ND  (g/sec) . 
mass  to  be  lofted  for 
i (g)  . 

total  mass  lofted  for  particle  size 
i for  all  time  (g) . 

Note  that  in  the  present  calculation 
aerodynamically  lofted,  and  only  one 
printed.  Other  calculations  using  four  particles  would  have  the 
values  for  each  particle  size,  printed  in  columns  below  the 
first  line  of  particle  1 values. 


only  one 
value  of 


particle  size  war 
this  variable  is 
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BEST  AVAItABLEXCrY 


Particle  Impact  Information 


OLD  CTCLl* 


ftCSTART  DUMP  AT 
OLD  CTCLC* 


OLD  CTCLL*  MSA< 


OLD  CTCLCa 


OLD  CTCLC«  IS3*< 


OLO  CTCLC*  IS02 


OLD  CYCLE*  1*12* 


OLO  CYCLE* 


OLO  CYCLE* 


OLO  CYCLE* 


OLO  TIHC* 

NO  or  INPAC 
NO  OF  PANTI 
PARTICLES  K 

T« 

OLO  T|HE« 

NO  or  IMPAL 
NO  OF  Parti 
PARTICLES  r 
OLO  TiHE* 

NO  OF  IHPaC 
NO  Of  PARTI 

Part iCLEs  A 
OLO  T|H£a 
NO  CF  INPAC 

NO  OF  Parti 

PART ICLES  K 
OLO  TIHE« 

NO  OF  INPAC 

NO  OF  Parti 
PARTICLES  K 
OLD  tine* 

NO  OF  INPAC 

NO  OF  Parti 
Particles  k 
OLO  T|HE* 

NO  OF  IHPAC 

NO  OF  Parti 
Particles  < 
OLO  TINE* 

NO  OF  INPAC 

NO  OF  Parti 
Particles  k 

OLD  TINE* 

NO  or  INPAC 
NO  or  PARTI 
PARTICLES  F 
OLO  TINE* 

NO  OF  INPAC 


S*2|t|*00  NCR  CYCLE* 
TS*  14 

CLES*  7 

ILLCO*  I 

00 

3.)V4a*00  NE«  CYCLE* 
TS*  It 

CLES*  7 

ILLLO*  0 

O.S7S4«00  NL»  CYCLE* 
TS*  20 

CLES*  S 

ILLCO*  I 

J.7StN*00  NCR  CYCLE* 
TS*  24 

CLES*  4 

ILLEU*  0 

H. 0224*00  NCR  CYCLE* 
TS*  26 

CLtS*  t 

ILLCO*  I 

N*2907*00  NCR  CYCLE* 
TS*  24 

CLES*  7 

ILLCO*  0 

N«N49S*00  NCR  CYCLE* 
rs*  27 

CLES*  7 

ILLCU*  0 

NCR  CYCLE* 
TS*  00 

CLES*  n 

ILLEO*  0 

S*00S9*00  NCR  CYCLE* 
TS*  07 

CLES*  II 

ILLED*  0 

S*040N*00  NCR  CYCLE* 
TS*  97 


NCR  CYCLE* 
24 

4 

0 

NCR  CYCLE* 

26 

9 

I 

NCR  CYCLE* 

24 

7 

0 

NE«  CYCLE* 
27 
7 
0 

NER  cycle* 
00 
1 1 
0 

NCR  CYCLE* 
07 
I I 
0 

NCR  CYCLE* 


1920.  NCR  tine*  0*0946*00  OT* 
TOTAL  NO  OF  INPaCTS* 
total  no.  of  PAKTICLCS* 
total  particles  killed* 

|9S4.  NCR  TINE*  O.S7S4*00  DT* 
total  no  of  IHPaCTS* 

TOTAL  NO.  OF  PARTICLES- 
total  particles  killed* 

1969.  NCR  TIME*  0*7S99*00  DT* 


TCT 

AL 

NO 

OF 

impacts* 

230 

TOT 

AL 

NO 

• OF 

PARTICLES* 

39 

1 

TOT 

AL 

PART  IC 

LES  KILLED* 

3 

6. 

NER 

1 T 

IHE* 

•:.0226.00  OT* 

1 

• ' 

TOT 

AL 

NO 

OF 

impacts* 

2S 

6 

TOT 

AL 

NO 

. OF 

particles* 

39 

7 

TOT 

AL 

PAI 

fIT  IL 

LES  KILLED* 

3 

1S62*  NEA  TIME*  H*2907*00  DT* 
total  no  OF  IHPaCTS* 
total  NO.  or  particles* 
total  particles  killed* 

1412.  NEA  TINE*  >t.i69S*00  DT* 
TOTAL  NO  or  INPaCIS* 
total  NO.  OF  particles* 
total  PAHTICLfS  KILLED* 

I4S9.  NE«  TINE*  9.7)77*00  OT* 
TOTAL  NO  or  INPACTS* 

total  no.  of  particles* 
total  particles  killed* 

1494.  NER  tine*  S.00S9.00  DT* 
total  no  of  impacts* 

TOTAL  NO.  OF  PARTICLES* 

total  particles  rilled* 

1796.  N£R  tine*  S.34)9.00  OT* 
total  no  OF  impacts* 

total  no.  of  particles* 
total  particles  killed* 

1609.  NCR  time*  S.8I09.00  DT* 
total  no  of  impacts* 


i.reTs^oi  shock  position* 

191 

)eo 

2 

1.7676-01  shock  POSITION* 
21U 
)67 
2 

1.7676-01  Shock  position* 


7676-01  shock  position* 


2.4617-01  SHOCK  POSITION* 
269 
90S 
9 

2*6617-01  SHOCK  POSITION* 
310 
9|2 
9 

1.7676-01  SHOCK  POSITION* 
337 
9|9 
9 

2.4617-01  SHOCK  POSITION* 
370 
930 

9 

2.6617-01  SHOCK  POSITION* 
907 
99  1 
*1 

3.S7S6-01  SHOCK  POSITION* 
9S9 


OLD  CYCLE  - cycle  number  of  previous  cycle  of  hydro  input. 

NEW  CYCLE  - cycle  number  of  subsequent  cycle  of  hydro  input. 

OLD  TIME  - time  of  old  cycle  (since  detonation) , sec. 

NEW  TIME  - time  of  new  cycle  (since  detonation),  sec. 

DT  - time  between  present  old  and  previous  old 

cycle,  sec. 

SHOCK  POSITION  - radial  location  of  air  shoc)c  at  old  time,  cm. 

NO.  OF  IMPACTS  - number  of  particle  impacts  in  dt. 

NO.  OF  PARTICLES  - number  of  particles  created  in  dt. 

PARTICLES  KILLED  - number  of  particles  deposited  in  dt. 

TOTAL  NO.  OF  IMPACTS  - total  number  of  particle  impacts  re- 
corded between  0 and  old  time. 

TOTAL  NO.  OF  PARTICLES  - total  number  of  particles  created  be- 
tween 0 and  old  time. 

TOTAL  PARTICLES  KILLED  - total  number  of  particles  deposited  be- 
tween 0 and  old  time. 
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5,  Sweepup  Dust  Cloud  Structure 

a? 
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• 
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0*0000 

0*0000 
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T 
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0*0000 
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2*4441«07 
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0*Oroo 

1 1 

4*5000*04 
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3*2441*02 
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0*0000 

12 
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3.4243*03 

0*C000 

0*0090 

O.onco 

13 
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5*2744*07 

5*2744*07 

o.ocoo 
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4*3650*03 

0*0000 
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0«OOP9 
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4*cr0i}*d4 

o.oooc 

c.ooco 
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0*0030 

O.OOOO 

IS 

4*5000*04 

1 *2346*04 

1 *2346*04 

0*0000 

o.ocoo 

o.onoo 

1 .0067*04 

0*0000 

0*0000 

0*0000 

14 

4*0000*04 

3*4927*07 

3*4927*07 

0*0000 

0*0C03 

0*0000 

5*7997*03 

0*0000 

0*0000 

0*0000 

IT 

7*5000*04 

4*42<0*C7 

4*4220*07 

0.0000 

0*0000 

0*0000 

1 *0003-04 

0*0000 

0*0090 

o.ooco 

It 

1*0003*05 

4.5454*07 

4*5659*07 

o.ocoo 

0*0000 

0*0000 

4.4439*02 
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0*0909 

0*0009 
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1 *0500*05 

1 *0044*04 

1 .COLA.Ct 

o.ooco 

0*0000 

0*0000 

1*0077*04 

0*0000 

0*9000 

o.cooo 

20 

1 • 1000*05 

2*62/|«(.4 

2*6221*04 

0*0000 

Q.COuO 

0*0000 

7*4969*03 

0*0000 

0*0030 

0 .o^*on 

21 

1 * 1 500*05 

1*21  70*04 

1 *2170*04 

0.0000 

0*0000 

0*0000 

3.5994*03 

0*0000 

o*ooco 

0.0''C9 

22 

1 * 2C00*05 

2*5215*04 

2*5215*041 

0.0000 

O.OUOO 

0*0000 

1 .7439*04 

0*0000 

0*0000 

0*0090 

23 

1 *2500*05 

2*2404*04 

2*2604*04 

0*0030 

0*0000 

0*0000 

1 *4532*04 

0*0000 

0*0090 

o.ooco 

24 

1 * 3CC3*05 

7*0147«04 

7 • 0 1 4 7*0tt 

0*0000 

0*0009 

0.0000 

1 *4529*04 

0*0000 

o*ooco 

o.ooco 

2S 

1 *3500*05 

2*3244*14 

2*3249*04 

0*0000 

0*0009 

0.0000 

1 • 1544*04 

o.cooo 

0*0090 

o.onoo 

24 

1 .4003*05 

4*  lv09*CA 

4* 10C9*0« 

0*0000 

o*ooco 

0*0000 

1 .1295*04 

0*0000 

0*0000 

o*nrc9 

2T 

1.4503*05 

4*4940*04 

4*6990*04 

0*0090 

c*ooco 

0*0000 

9 •* 1 26*03 

0*0000 

0*0000 

o.onoo 

2t 

1 *5CC0«05 

1 *0174*04 

1 *0179*09 

0*0000 

o*ooco 

0*0000 

4*3042*03 

0*0000 

o.cooo 

0*0090 

24 

1*5500*05 

4*7  1 4h*04 

4*7  194*06 

0*0000 

0*0000 

0*0000 

7*4149*03 

0*0090 

0*0090 

o.ooou 

30 

1 .4000*05 

1 • 1 742*04 

1 * 1 792*09 

o.ooco 

o.ooco 

0*0000 

7*5549*03 

0.0000 

0*0030 

0.CO90 

31 

1 .4500*05 

4*7454*04 

4*7454*011 

0*0003 

0.0000 

0*0000 

4*2439*03 

0.0000 

0*0090 

Ci.QOOO 

32 

1 * 7000*05 

1 *0234*04 

1 *0239*09 

0*0090 

0*0000 

0*0000 

5*  3443*C3 

0.0000 

0*0090 

0*000.9 

33 

1 *7500*05 

4.0295*04 

9.4295*04 

9 *0000 

0*0000 

0.0000 

4.7917  *03 

0*0000 

0*0090 

o.ooco 

34 

1 *4C00*05 

7 • 1 3*2  *U4 

7 * 1 342*04 

o.ooco 

0*0000 

0*0000 

4.5437*03 

0*0090 

0*0000 

C • 0000 

3S 

1 .4500*05 

7*5443*04 

7.5463*04 

0.0000 

0*0000 

0.0000 

4*0423*03 

0.0000 

0*0990 

0*0000 

34 

1 .4000*05 

2*3223*09 

9*2229*06 

0*0000 

0*0000 

0.0000 

3.7944*03 

0*0000 

0*0900 

O.OOOO 

37 

1 *4500*05 

4*3471*09 

4 • 6 7 1 2*04 

0*0000 

0*0003 

0.0000 

3*3003*03 

0.0000 

0*0009 

0*0003 

34 

2*0000*05 

2* 3544*09 

4*3441*04 

0.0000 

0*0000 

0.0000 

3*2310*03 

0.0000 

o.ooco 

o.onoo 

34 

2.Q5oO*05 

1*7341*09 

1*0941 *09 

0*0000 

o.ooco 

0.0000 

3.032 1 *03 

o.cooo 

o.ooco 

0*0000 

40 

2*  1000*05 

5*7975*04 

5.7975*04 

o.ocoo 

0*0900 

0*0000 

2.5440*03 

0*0009 

o.oooc 

O.OOPO 

4 1 

2* 1500*05 

7*0259*04 

7*0259*08 

0.CC90 

0*0000 

0*0000 

5*3*05*02 

o.onoo 

0*0090 

0*0000 

42 

2*2000*05 

5*4754*04 

5*4754*04 

0*0000 

0.0000 

0.0000 

2*2442*03 

0*0090 

0*0000 

0 * 3000 

43 

2.2500*05 

7.5051*04 

4*7051*04 

o.cooo 

0*0000 

0*0000 

1 *9005*03 

0.0000 

0*09C0 

0.0000 

44 

2*3000*05 

4*0006*04 

4*0004*04 

0.0000 

0*0003 

0.0000 

1 .6303*03 

0*0090 

o.ooco 

0*0000 

4S 

2*3500*05 

0*0L0C 

0*0000 

0*0000 

0.0000 

0*0000 

0*9000 

0*0000 

0.0900 

o.ooco 

44 

2.4C00*05 

2*5444*06 

2 *5444*04 

0.0000 

0*0090 

0*9000 

1*2121*03 

0.0000 

0*0030 

0.0O09 

47 

2*4500*05 

3 . 2044  *04 

3*2044*04 

o.cooo 

o.ooco 

0*0000 

9*CA99*02 

0*0000 

0.0000 

o.ooco 

44 

2*5000*05 

4 .4474*04 

4.4474*04 

o.ocoo 

0*0000 

0*0030 

2*  74 1 7*02 

9*0000 

0*0900 

o.cocu 

44 

2*5503»C5 

O.OCOQ 

0*0900 

0*0000 

0*0003 

O.OQOO 

0*0000 

0.0000 

o.oooc 

o.uroo 

SC 

2.4000*05 

0.C60C 

O.CLOO 

o.ccoo 

0*OoOO 

0*0000 

0*0000 

0.0000 

0.0000 

o«oroo 

S 1 

2*4500*05 

0.0 uno 

0.C3CC 

o.cooo 

0*0000 

C *0000 

0.0000 

0*0000 

0*0000 
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J 

T 19 

tOlMAS 

HA55  neOTt 

SUMS  1 

S04S2 

SUMS! 

SUNS4 

1 

5*0000*01 

3* 1 344*09 

2*2400*10 

3*  1 366*09 

0*0000 

0.0000 

0*0000 

2 

1 *OCC0*O2 

7.4912*04 

2*  1 731  ^iu 

7 * 4f i 2«0A 

0*0000 

0.0000 

0*0000 

3 

2*0300*02 

4.4  9<t6*0A 

2*0447*  Ij 

4*4946«U4 

0*0000 

0.0000 

o.ocoo 

4 

3*0033*02 

7.C934*0A 

2*0172*10 

7*0936.04 

0*0300 

0.0000 

0.0000 

s 

4.0003*0/ 

5*6467*09 

1 *4404*10 

5*6467«04 

0*0000 
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4 
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0*00UC 

1 .9404*10 
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7 
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0*0300 
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4 
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4 
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10 
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0*0000 

1 2 
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0*030C 
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0*0000 
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0*0000 

•2*0400*04 
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For  time  T,  total  mass  of  aerodynamically  lofted  dust 
in  the  grid  is  printed,  with  mass  of  dust  for  each  individual 
particle  size  also  indicated. 


PKILLl,  PKILL2 
RZ  (I) 

TOTMASS (I) 
SUMSl (I) 

SUMS2 (I) 

SUMS3 (I) 

SUMS4 (I) 

YMAXl (I) 

YriAX2  (I) 

YMAX3 (I) 

YMAX4 (I) 

YZP (J) 

TOTMASS (J) 

MASS  ABOVE 
SUMSl (J) 

SUMS2 (J) 

SUMS3 (J) 

SUMS4 (J) 


- amount  of  mass  deposited  (using  two  particle 
deposition  routines) . 

- X coordinate  of  center  of  zone  i. 

- total  dust  mass  in  zone  i for  all  j. 

- mass  of  dust  particles,  size  1,  in  zone  i 
for  all  j. 

- mass  of  dust  particles,  size  2,  in  zone  i 
for  all  j. 

- mass  of  dust  particles,  size  3,  in  zone  i 
for  all  j . 

- mass  of  dust  particles,  size  4,  in  zone  i 
for  all  j . 

- y coordinate  of  maximum  height  reached  by 
particle  size  1. 

- y coordinate  of  maximum  height  reached  by 
particle  size  2. 

- y coordinate  of  maximum  height  reached  by 
particle  size  3. 

- y coordinate  of  maximum  height  reached  by 
particle  size  4. 

- y coordinate  of  top  of  vertical  zone  j . 

- total  dust  mass  in  zone  j,  for  all  i. 

- total  dust  mass  above  zone  j for  all  i. 

- mass  of  dust  particles,  size  1,  in  zone  j 
for  all  i. 

- mass  of  dust  particles,  size  2,  in  zone  j 
for  all  i. 

- mass  of  dust  particles,  size  3,  in  zone  j 
for  all  i. 

- mass  of  dust  particles,  size  4,  in  zone  j 
for  all  i. 
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6.  Grid  Dust  Concentrations  (Aerodynatnically  Lofted  Dust) 
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7.  Concentration  Level  Printer  Plots  — All  Dust 
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8.  Concentration  Levels  for  Particle  Size  1 
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9.  Dynamic  Pressure  Levels 
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10.  Horizontal  Dynamic  Pressures 
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11.  Preloaded  Particles 

The  entire  retinue  of  printout  formats  are  available 
for  the  preload  particles  only.  The  following  two  plots  are 
provided  for  the  concentration  and  horizontal  dynamic  pres- 
sure of  1 mm  diameter  crater  particles  at  t = 3.2181  secs  in 
a 1 MT  surface  burst. 


APPLE  27 


concentration  levels. 

1 . saoo«06 

I .0030*06 

s.oouo*os 

3.aooQ*os 

2. Q00Q*QS 
1 ,0000*05 
5.0000*0“* 

3.0000*0“* 

2.0000*0“* 

I . 0000*0“* 

5.0000*03 
3.0000*03 
2.0000*03 
I .0000*03 
5.0000*02 
“*.0000*02 
3.0000*02 
2.0000*02 
I .0000*02 
5.0000*01 


I. 00-03  T*  3.2181*00  HALF  OECAOE  INTEGER  PLOTS 


3321  I 
3“*  “*3 


I . I 000*05 


♦ 

2.1000*05 


3.1000*05 


APPLE  27 


HORiz  dtnam  pres  e« 

I .5000*06 
I .0000*06 
5.0000*05 
3.0000*05 
7.0000*05 
I .0000*05 
5.0000*0“* 

3.0000*0“* 

2.0000*0“* 

1 .0000*0“* 

5.0000*03 
3.0000*03 
2.0000*03 
I .0000*03 
5.0000*02 
■*.0000*02 
3.0000*02 
2.0000*02 
I .0000*02 
5.0000*01 


I .00*05  T>  3.2181*00 


I . IOOO*OS 


FULL  OECAOE  INTEGER  PLOTS 


6655  5 

6666 


♦ ♦ 

2.1000*05 


3.1000*05 


136 


APPENDIX  B 


BLAST  WAVE  FLOW  FIELDS  FOR  SCOUR2  CALCULATIONS 

MIDDLE  GUST 

To  compute  the  dust  cloud  development  behind  the  blast 
wave  from  the  100  ton  MIDDLE  GUST  detonation,  the  flow  condi- 
tions at  the  outer  edge  of  the  boundary  layer  were  computed 
on  the  basis  of  an  idealized  flow  model.  An  analytic  fit  to 
the  air  shock  position  as  a function  of  time  (X^(t))  and 
values  for  the  duration  of  the  positive  phase  at  a given 
radial  location,  At_^(X^),  are  the  basic  determinants  of  the 
flow  field.  For  MIDDLE  GUST,  a fit  to  the  data  in  Table  B.l 
from  Reference  B.l  yields 


x^Ct) 


= 170 


/t 


h) 


0.525 


(B.l) 


Xs(t) 


= 4.46 


i I 

I 20  ( 


-475 


(B.2) 


Values  for  At_^  are  taken  directly  from  Table  B.l. 

Since  shock  location  and  speed  are  known,  we  may  com- 
pute the  air  velocity,  density,  and  internal  energy  behind  the 
shock  from  the  shock  relations  for  an  ideal  gas  with  specific 
heat  ratio  of  1.4.  This  gives  us  the  peak  values  of  these 
parameters: 


and  (E  , ) , where  N 
peak 

designates  the  zone  number  where  the  shock  is  located;  i.e.. 


^^peak^j^ 


^^peak^„ 


X - ^ < 
^N  2 


^S  " 


X + ^ 
^N  2 


(B.3) 


where  AX  is  the  zone  size  (2.5  m for  MIDDLE  GUST).  The  flow 
parameters  for  the  100  ton  MIDDLE  GUST  blast  wave  are  pre- 
sented in  Table  B.2. 
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Table  B.l.  Predicted  Air  Blast  Parameters 


for 

MIDDLE 

GUST  Event 

III 

Ground 

Range 

(ft) 

Arrival 

Time 

(msec) 

Shock 

Pressure 

(psi) 

Positive 

Duration 

(msec) 

Positive 

Pressure 

Impulse 

(psi-^nsec) 

Dynamic 

Pressure 

(psi) 

Dynamic 

Pressure 

Inpulse 

(psi^nsec) 

17 

0.80 

20,000 

9.0 

48,000 

130,000 

150,000 

18 

0.85 

15,000 

9.0 

35,000 

110,000 

130,000 

20 

0.94 

10,000 

8.9 

21,000 

78,000 

110,000 

22 

1.0 

5,000 

8.8 

14,000 

60,000 

91,000 

25 

1.2 

3,000 

8.8 

8,000 

43,000 

74,000 

30 

1.5 

2,000 

8.8 

4,000 

18,000 

57,000 

70 

4.7 

1,000 

8.3 

1,000 

8,000 

19,000 

92 

7.1 

750 

8.1 

970 

4,700 

14,000 

120 

11.0 

500 

7.9 

1,000 

2,800 

10,000 

170 

20.0 

250 

22.0 

2,200 

1,100 

9,500 

250 

40.0 

100 

62.0 

1,200 

310 

5,100 

320 

60.0 

50 

81.0 

830 

75 

2,700 

380 

90.0 

30 

95.0 

640 

18 

1,200 

440 

115.0 

20 

110.0 

560 

8 

350 

560 

195.0 

10 

130.0 

440 

2 

72 
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Table  B.2.  Idealized  Flow  for  MIDDLE  GUST  II  Dust  Lofting 
Calculations  Flow  Parameters  Just  Aft  of 
the  Air  Shock  as  the  Shock  Enters  the 
Various  Zones  — MIDDLE  GUST  II 


r: 

N 

TShOCK 

VSnOcK 

Snockn 

^SnOcK 

^SmOCK 

C5W0C* 

OShOCI 

U^CaC 

OCCTS 

:0\E  LO- 

ZONE 

ARRIVAL 

SHOCK 

SHOCK 

SHOCK  over- 

SHOCK 

SHOCK 

SHOCK  TEM- 

PEAK AIR 

POSITIvl 

CATION,.'! 

NO. 

Tl.M£,SEC 

SPLhD, 

HACH 

pressure 

DENSITY 

ENERGY, 

PERATURE, 

velocity. 

PHAEE  DURA 

CM/SEC 

NUMBER 

(ERCS/CC) 

(CM/CC) 

(ERCS/CC) 

•k 

(ovsr.c) 

TION.SLCS. 

3*1414»o3 

3«2| « 1 *05 

4.6757 

00 

4.O0p0*O7 

6,0004-03 

3*7320* 10 

5,2650*03 

2.6527*05 

0*orc2 

:i.:s 

2 

4.otl4-o4 

2*e420«oS 

•*6437 

00 

7*3350*07 

5*43*1-03 

3*05| 3* 10 

4. 304  7*0  ) 

2*3704*05 

1*0200-32 

:j.7s 

> 

4,4433-03 

1*4246*05 

7.657a 

00 

6.0541*07 

5.0501-03 

2*5536*10 

3.6024.03 

2*  1560*05 

♦ .5000-0 

2o . 2 S 

0 

S,4t37-03 

1.4 1 2S*o5 

7*2006 

00 

5*0473*07 

5,7771-03 

2*1766*10 

3.07)6*03 

1.4717*05 

9,  1200-03 

28.  ‘S 

6 

7. 0413-03 

l*2244*eS 

6*6624 

CO 

4*3526*07 

5.6422-03 

1*6602*10 

2.66  36*03 

1,0161*05 

• . ICOO-O) 

31 . 2S 

4 

••12SS-03 

1*0741 *os 

6,1471 

00 

3*7636*07 

5.6030-03 

1*6565*10 

2.3)’0*O3 

1 .6e)i*C5 

0.2500-03 

35. :s 

47^8-03 

1 .434**cS 

5,7455 

CO 

3.7647*07 

5. 5127-03 

1*47)6*10 

2.0704*03 

1 ,56^2*05 

7 , 4C2C-0) 

58. :s 

• 

» •CSC3*C2 

1 .6222-05 

5 , H14  6 

CO 

2.4020*07 

5.4143-03 

1 • 322  )*  1 0 

1 *6655*03 

i.467j*05 

7 , 4Ci«- S J 

38. '5 

♦ 

1 *21  1 4-C2 

1 • 7 1 ea -cs 

5*1354 

00 

2*5606*07 

5.3241-03 

1 * 1 464*  1 0 

1 .6f 15*0) 

1 •)  7 a 1 *05 

6 . 5 2 2 0 • 2 ) 

41. 2S 

|0 

1*371 )-Cl 

1 . 42  7 |«0S 

4,0616 

00 

2*3|  10*07 

5.2275-03 

1 *04  16*10 

1 .5101 *0) 

1 .2966*05 

1 , 223C-C2 

43.75 

1 

1 tsirt-ci 

1 .ShS 1 *05 

4,6166 

CO 

2*0625*07 

5. 1 799-03 

1 .0024*  to 

1.4  1 42*0) 

1 • 22  7 2*  05 

1 .7930-02 

46.25 

I 2 

1 ,44  4 7-02 

1 .4  7 1 3*0S 

4,34*1 

00 

1 . 067  1 *o7 

5.031 7-03 

4 . 26 1 2*04 

i . 3065*03 

1.(627*05 

2.  J6  ';&-02 

4S . :s 

1 3 

I *4664-02 

1 * 404**0S 

4,1470 

00 

1.7166*07 

4.4331-03 

6*6027*04 

1.21)7*0) 

1*1041*05 

2.v2:o-02 

S 1 . 2S 

1 4 

2.CSC4-02 

1 . J440-0S 

4.0156 

00 

1 .5  7 22*07 

4,e3i*-oj 

t*0J05*n4 

1 . I 324*03 

1 , 0 5 C 5 * 0 5 

3. 4930-22 

S3.7S 

IS 

2 ,2404-02 

1 • 26  6&-CS 

3*6503 

CO 

1 .4111*07 

4,73*3-03 

7 ,5297*04 

1*0623*03 

1*0011-05 

4.P6C0-C2 

So.  2$ 

1 * 

2.43t4-02 

1 • 2 3 70*OS 

1,6405 

CO 

1 • 3 3 1 4 *07 

4 , 6 )«0-03 

7*0607*04 

1 .C03i*0i 

4 « 56  07  *04 

4,6)C«1-02 

SS.'S 

\ 7 

2,444S-02 

1.1410*05 

3,5566 

00 

1 .2317*07 

4*5407-03 

6,6463*04 

4,1196*02 

4 , 1 1 1 1 • C4 

5.2'‘.4L-02 

61 . 25 

1 • 

2,aS6)-02 

1 « I 4 7* • QS 

3,4247 

00 

1 . Ii30*07 

4.4112-03 

6 , 350  7 *07 

0.45’5*O2 

6*  752  j *Ci 

5.)9C0-CJ 

65.75 

t « 

3,ceco-02 

1 * 1 0 7 6 -qS 

3,3102 

00 

1 .0637*07 

4 , )4«d-03 

6,0400*0’ 

6.5211*02 

6 . )0’ 2 • 04 

5.573C-C2 

06. 2$ 

20 

3*  30*4-02 

1 . 0706*05 

3*  1490 

00 

♦.4256*06 

4,2546-03 

5*  ?6Ca*0’ 

« * 1 2 72*02 

t.05SC*04 

5. 760L-02 

63.:’S 

7 1 

3.S470-02 

1 .0360*05 

3,0455 

CO 

4.2617*06 

4.1617-03 

5,5q40*C4 

7 , 7 7 20*  02 

7 . 7 3 2 1 *C4 

5.9530-02 

71.25 

12 

3.74  2 1 -02 

1 .0036*05 

2,4467 

00 

6,7046*06 

4.0702-03 

5.2610*0’ 

7 ,4504*02 

7.13)2*04 

6 . 1 IOC- 02 

■3.'S 

73 

4.04SI-P2 

4.?  320*o'' 

2.4062 

00 

6 , 1 765*06 

3. 4*02-03 

5,a/)’*0’ 

7. 1502*02 

7.1516*04 

4.)?CC-02 

‘6.25 

2 4 

4.3CSt-02 

••4.,fi3*04 

2,6232 

00 

7.6942*06 

3-0410*03 

4,6851 *0’ 

6 « 6 9 1 6 *C2 

6«ee57*Ci 

4.5l0:-02 

‘8.75 

IS 

4.S742-02 

4 . 1 607*04 

2,7i32 

00 

7.2615*06 

3.eo5l-C3 

4 . 7 121*0’ 

6.616I  »02 

6 • 6 ) 39  *04 

6 . 7 COO-02 

SI . 25 

/4 

••.*503-02 

• .4;b?*o'» 

2 . 6*  7 1 

00 

6*0606*06 

3.7201*03 

4 ,55h0  *0’ 

6,4247*02 

6.3952*04 

6 . 7 700-  02 

83.75 

27 

4,  1 34  1 -02 

6,6407*04 

2.54*6 

00 

6.49^5*06 

3.6)60-03 

4.4003*0’ 

6.2l42*C2 

6 . 1 60  )*C4 

6.«SC'j-02 

50.25 

2$ 

S,  41S4-02 

*,4*S6«04 

2.524* 

DO 

6*  15)5*06 

3.555)-0) 

4,2710*0’ 

4 . 02  96  *02 

5,9V2)*51 

6,42Pu-02 

88.  75 

S » 7 2 4 7 -C2 

0.2527*04 

2.4*54 

00 

5»6h07*o6 

3.4755-0) 

4 • 1 49  8 *0’ 

5.05*1*02 

S.7i6)*C4 

7.0''?  -02 

91.25 

30 

4. 0314-02 

0 . 0506*04 

2.4055 

00 

6*5511*06 

3.)4;6-0) 

4.0)47*04 

5.6920*02 

5, 5*95.04 

7.070O-C2 

Q . 7 5 

3 t 

4.34S7-C2 

7.0587*04 

2.3462 

00 

5*26)2*06 

3.)2l5-03 

3,4270*0’ 

5.51  1 2*02 

5 . )6 1 2 * C4 

7 , 1 5 J C • 0 2 

96.25 

)2 

4.447S-C2 

7.6762*04 

2.2437 

00 

5.0342*04 

3.2472-03 

3*0263*0’ 

5.4039*02 

5 • 1 0C9  *04 

>,22CO-02 

9S.‘5 

33 

4,44*4-02 

7 « Sc  2 3 * 0 4 

2.2417 

00 

4.6025*06 

3.  t 716-03 

3 . 7)55*  0’ 

5.270C*02 

5.0276.04 

7 . j:oc-02 

101.25 

3 ** 

7 , 3)34-c2 

7.33*6*04 

2.1422 

00 

4.5665*06 

3 . 1 039  “C) 

3*6100*0’ 

5.1177*02 

1 . 8 * 1 6 • 0" 

7,4233*02 

1 0 3 . ‘ 5 

3S 

7 , 474.I-C2 

7.1764-04 

2*1411 

00 

4,3016*06 

3.0)49-03 

3*5676*0’ 

5.0)31*02 

4.66(7*04 

7.55CC-02 

lOo. 25 

3* 

•,C3C3-02 

7*0272*04 

2.0447 

CO 

4,1962*06 

2,4677-03 

3.1415*0’ 

4,9257*02 

4,5277*04 

7,67CC-02 

108.75 

3T 

t,  3*40-02 

6.8625*C4 

1.05&5 

00 

4*0145*06 

2.4022*03 

3.4199*0’ 

4,62*6-02 

4. 3793*04 

7,  0CC':-C2 

111.25 

)• 

•.7S47-02 

4 * 7 4 4C  *04 

2,0l5| 

00 

3*0537*06 

2*0361-03 

3,3526*0’ 

4.7293*02 

4.2)60*24 

7.920L-02 

115.75 

4.131  1-02 

4.411 2*04 

1.4754 

PO 

3,6474*06 

2 . 7 76  2-03 

3.2092*0’ 

4 , 6103*02 

1.0976*01 

6.05Cr-02 

116.25 

• 0 

4.SI24-02 

4,4036*04 

1,4374 

00 

3*5514*06 

2.7157-03 

3.229)*0’ 

4,5556*02 

3.96)4*04 

6.  1 7C0-02 

113.75 

1 

4,4022-02 

4*  34  1 5*04 

1 *4006 

00 

3,4134-06 

2.6566*03 

3. 1 724*0’ 

4, 4 759*02 

3.6)10*04 

6.3CCC-C2 

1 : 1 . : 5 

« 2 

1 »C244-C 1 

4,1434*04 

1.0657 

CO 

3.2633*06 

2.5491-03 

3, 1 l9p.o’ 

4.4002*02 

3»7C0i*C4 

6 . 4 1 CC-02 

123.75 

43 

1 *0703-01 

4.  1 300*04 

1*0314 

00 

3,1604-06 

2. 54)6-03 

3*0661 *0’ 

4,3205*02 

3.5666*04 

•.52:0-02 

126.25 

4 <* 

1 * t n 4-01 

6.021^*04 

1 , 794  1 

00 

3,0412*06 

2.4042-03 

3*0l”*0’ 

4,2601*02 

3.4603*04 

6.61C0-C2 

i2o.  5 

44 

1 . |Si  j-Ci 

S . 7 1 70  *04 

1 , 766(1 

oa 

2.4342*06 

2.  4)63-03 

2,4740»0’ 

4 1 1 956  *02 

) . JS  )5*C4 

6 , 750-'-O2 

151.25 

44 

1 . 1 4S4-0 1 

S . 6 1 5 4 *04 

1*7370 

00 

2 . 0 300*06 

2.  3649-03 

2*’)0)*0’ 

4 , 1 340*02 

).  2 1 1 1 *24 

8,6602-02 

I 5 3 . ‘ 5 

42 

1 . 234  1-0 J 

S . 7 j 6 4 *0  4 

I.7CC7 

CO 

2-73  1 7*96 

2.3346-03 

2.6H66*0’ 

4.0752*02 

) . 1 ) ) 1 *Ci 

6 . ’ 7 CO-0  2 

156.25 

4l 

i .2t3)-0l 

5,4243*04 

1 ,6605 

00 

2*4375*06 

2.2e( 1-03 

2.8109*09 

4.0192*02 

).027)*C4 

9.0100-02 

158. ‘5 

4» 

i . 32(i2'Oi 

5,5333*04 

1*6531 

00 

2*5463*0* 

2.2367-03 

2.6lC’*0’ 

3.4656*02 

2.424)*04 

9.2:00-02 

141.25 

SO 

1 * 3737-01 

5,4654*01 

l«627| 

00 

2*4636*06 

2.14.5-03 

2*7716*0’ 

3.4141*02 

2*02)6*04 

9.3000-02 

143.75 

SI 

! *4200-01 

S*340«*0« 

1*6017 

00 

2-3024«C6 

2. 1476-03 

2.7344*04 

),I6S4*02 

2,7256*04 

9,4300-02 

146. 25 

SI 

1 *4*70-01 

5*1761 *04 

1*5771 

00 

2*306|«06 

2*1034-03 

2*706S*0’ 

3*6  103*02 

2*6301*04 

9*5000*02 

148. '5 

13 

1 .S  14  7-C  1 

5, 1465*04 

t .5533 

00 

2*1326*06 

2.0611-03 

2,6744*04 

3*7732*02 

2*5)66*04 

9.6000-02 

151.25 

S4 

1 •S431-0I 

5*111  3*04 

1 *5303 

00 

2*  1626*0* 

2.0201-03 

2*6434*0’ 

3. 7244*02 

2*4452*04 

’.70CO-02 

153. ‘5 

ss 

1*4)  2 3-CI 

5*0465*04 

1 *5o74 

00 

2*0460*06 

1 *4746-03 

2*6|4)*0’ 

)*6002*Q2 

2.3559*04 

9,6000-02 

156.25 

S4 

1*442?-Q| 

4,474C*04 

1*4662 

00 

1*0322*06 

1 *4406-03 

2.565’*0’ 

3, 616  1 *02 

2. 268S*0« 

♦.♦OCO-02 

153.75 

si 

1 .7120-01 

4,4036*04 

1 * 4652 

00 

1*471 2*06 

1 .402S-C3 

2,55B5*C’ 

3.6045*02 

2 . 10  29  *04 

1 .COCO-Ol 

161.25 

St 

1 * 744  1 -01 

4,6353*04 

1 ,4116 

00 

1*4120*06 

1 . 0 65  3-0  3 

2,5)21 *0’ 

3,5722*02 

2.09’1*04 

1 ,C ( 20-01 

Io5. 75 

1.81*2*01 

4.7*40*04 

1 * 4250 

00 

1 •65«6«06 

1 . 6242*03 

2.5066*0’ 

3, 5)6). 02 

1,01  70*04 

1 ,C20C-C1 

166.25 

*c 

1 .0440-01 

4,7045*04 

1 • 405  7 

00 

1.0032*06 

1.7410-03 

2.4620*0’ 

3,5015*02 

1 .9)65.04 

1 ,C JCl-Cl 

loS.  *5 

* 

1 . 422S-0I 

4 , * I 1 4 • 04 

1 .3670 

00 

1.751  6*0* 

1 .7597  -03 

2.4562*0’ 

3.4679*02 

1.0575*04 

1 .OiCO-01 

1*1,25 

4l 

1 .47*7-01 

4.5610*04 

1,3606 

00 

1 • 7Q2h*06 

1*7261-03 

2*  4)5 1 *0’ 

3.4)54*02 

1 . ?0CO*C4 

J .0500-01 

1*3. *5 

43 

2.0314-01 

4.521 7*04 

1.3511 

CO 

1 *655  1 *06 

1 ,64)6*0) 

2.4126*0’ 

3.10)4.02 

1,70)4*04 

1 .OAO&'O  1 

i‘6.:5 

4 4 

2.0072-01 

4.4*41-04 

1.13)4 

00 

1 .6096*0* 

1 ,6622-03 

2. )9| 2*0’ 

3 . ) 7 )i*02 

1 .6292*04 

1 .0700-01 

1*8. ‘5 

4S 

1 . 1 4 3*-c  1 

4,4074*04 

1.3171 

00 

1 .565  0*06 

1,6)14-0) 

2* J7C2*3’ 

3 . ) 4 3 6 • 02 

1 .5556*04 

1 .0030-01 

131.25 

44 

2.2007-01 

4 • 35  3 3*04 

1 *3000 

00 

1 * 5230*06 

1 .60  1 3-0) 

2.3190*0’ 

3. 3 1 50*C2 

1.40)7*04 

1 -0000-01 

185. *5 

47 

2.2504-01 

4 . 3000*04 

1.2614 

00 

1 • 40  3 )*06 

1 ,5720-0) 

2* ))C0*0’ 

3 . 2 C-  7 1 * 0 2 

1.4127*04 

1 .C900-0 1 

150.25 

*t 

2*  3 144-01 

4*2461 *04 

1 * 2643 

00 

l.iiij^o* 

1 ,51)5-03 

2, 3 1 C 7*0’ 

3.2544.02 

1 . )i 29  * 04 

1 ,C4O0-CI 

IS8.  ‘5 

4* 

2*3741-01 

4.  1 4 75*04 

1.2542 

00 

1 *4060*06 

1*5157-0) 

2 . 2’ 1 4*0’ 

3.2))4.02 

l.274)*C4 

1 ,C9CO*0l 

191.25 

70 

1. 4360-01 

4,1462*04 

1 .2345 

00 

1 * 37C**C6 

1 *4606-03 

2*27)7*04 

3 • 20  76  * 02 

1 , 206  7 *04 

1 .0900-01 

1.45. *5 

71 

2. 44**. 01 

4 • t 000*04 

1.2251 

00 

1.3356*06 

1 .462  2*03 

2*2556*04 

3i  1 0 25*  02 

1.1402*04 

1 ,0400-01 

196.25 

71 

2iSS8C-0l 

4,0530*04 

i • 21  tn 

00 

1 « 302  1 *06 

1 *4364-03 

2.2)05*0’ 

3.1500*02 

1.0716*04 

1 • 1 coo-3  1 

194.  *5 

7 

2. *2:0-01 

4.0072*34 

1.1473 

00 

1 * 2694  *0* 

1*411  3*03 

2.2215*0’ 

3.1341*02 

1 *0100*04 

1 , 1 000-01 
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The  wind  field 
mated  in  the  positive 


at  a particular  location  may  be  approxi- 

[ B 2 1 

phase  by  the  following  prescription: 


where 


= peak  air  velocity  in  zone  N 


(B.4) 


(B.5) 


(B.6) 


shock  arrival  time  at 
zone  N 


At  ) = positive  phase  duration  in  zone  N 

'n 


X„  = radial  location  of  Nth  zone 

N 


E = internal  energy  of  air 

p = density 

e subscript  refers  to  conditions  at  outer  edge  of  the  boundary 
layer.  Spherical  divergence  in  the  flow  field  is  taken  into 
account  by 


u 


(B.7) 


The  negative  phase  is  prescribed  to  be  five  times  as 
long  as  The  horizontal  wind  velocity  is  constant  at 

0.1  and  the  air  density  and  energy  are  held  fixed  to 

the  values  computed  at  the  conclusion  of  the  positive  phase. 


MIXED  COMPANY 

The  hydrodynamic  flow  field  used  in  the  SCOUR  calcula- 
tions to  simulate  the  MIXED  COMPANY  (500  tons  TNT)  blast  wave 
flow  field  is  a scaled  version  of  the  MIDDLE  GUST  flow  field 
described  above.  The  location  of  a given  overpressure  loca- 
tion, the  time  of  air  shock  arrival,  and  positive  phase  dura- 
tion were  each  scaled  by  a factor  of  1.71,  the  cube  root  of  5, 
the  ratio  of  the  MIXED  COMPANY  yield  to  that  of  MIDDLE  GUST 
(100  tons  TNT) . A tabular  listing  of  flow  field  variables  is 
provided  in  Table  B.3. 
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Table  B.3,  Tabular  Listing  of  Flow  Field  Variables  — MIXED  COMPANY 
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Tech.  Info.  Center/S-1930 
R.  K.  Plebuch,  Rl-2078 
I.  E.  Alber,  Rl-1008 
Peter  K.  Dai,  Rl/2170 


TRW  Defense  8 Space  Sys.  Group 
San  Bernardino  Operations 

ATTN:  E.  Y.  Wong,  527/712 

Universal  Analytics,  Inc. 

ATTN:  E.  I.  Field 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

URS  Research  Coi.pany 

ATTN:  Technical  Library 

The  Eric  H.  Wang 

Civil  Engineering  Rsch.  Fac. 

ATTN:  Larry  Bickle 
ATTN:  Neal  Baum 

Washington  State  University 

Administrative  Organization 

ATTN:  Arthur  Miles  Hohorf  for  George  Duval 

Weidlinger  Assoc.  Consulting  Engineers 
ATTN:  J.  W.  Wright 
ATTN:  Melvin  L.  Baron 

Weidlinger  Assoc.  Consulting  Engineers 
ATTN:  J.  Isenberg 

Westinghouse  Electric  Corp. 

Marine  Division 

ATTN:  W.  A.  Volz 


